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Abstract. The existence of dark matter provides strong evidence for physics beyond the Standard Model. 
Extending the Standard Model with the Peccei-Quinn symmetry and/or supersymmetry, compelling dark 
matter candidates appear. For the axion, the neutralino, the gravitino, and the axino, I review primordial 
production mechanisms, cosmological and astrophysical constraints, experimental searches, and prospects 
for experimental identification. 

PACS. 95.35.+d Dark matter - 14.80. Mz Axions - 14.80.Ly Supersymmetric partners of known particles 



1 Introduction 



(Ci^B) whose detection is a major challenge; cf. [B] and 

references therein. While X?™"' ° ~ 0.003% for three 
massless neutrino species, we know from oscillation ex- 
periments that at least two neutrinos have small but 
non-zero masses, J2i ''^vi ^ 0.05 eV, which are not part 
of the Standard Model, such that > 0.1%. More- 
over, since light neutrinos are hot dark matter, which 
leaves an imprint on large scale structure (LSS) [J, one 
can extract the cosmological limit J^i ^ ^'(1 eV) 
or equivalently < 1.5% [8]. 
provided in the form of dark energy, non-baryonic dark (iv) The dominant amount of the non-baryonic dark matter 



Based on numerous cosmological and astrophysical stud- 
ies, we believe today that our Universe is flat and thus that 
the total energy density is critical, ptot — Pc = SHq/SttG^ 
or Htot = Ptot/Pc — 100%, with contributions of [T]|2] 



Ha 



a 



72%, i7din 

. ~ 0.005%, 



- 23%, c 

0.1% < < 



i 4.6%, 
1.5% 



(1) 



matter, baryons, photons, and neutrinos, respectively. Here 
fli = Pi/ Pc, Hq is the present Hubble expansion rate, and 
Gn is Newton's constant. The density parameter fi^y is 
given basically by the photons of the cosmic microwave 
background (CMB). The understanding of the remaining 
f^tot — — 99.995% requires physics beyond the Stan- 
dard Model: 

(i) The amount of dark energy f2A could be provided by 
a cosmological constant. However, the fact that the 
cosmological constant or vacuum energy inferred from 
a naive quantum field theoretical estimate exceeds Qa 
by 120 orders of magnitude is the most serious fine- 
tuning problem in physics. An alternative explanation 
of dark energy is a slowly evolving scalar field [SUUIl] 
that cannot be part of the Standard Model. 

(ii) The amount of baryons J7b — inferred from the CMB 
anisotropics, the framework of primordial nucleosyn- 
thesis, and the observationally inferred abundances of 
light nuclei — indicates a matter- antimatter asymmetry 
that cannot be explained within the Standard Model. 

(iii) The amount of neutrinos fi^, is understood to be pro- 
vided in the form of the cosmic neutrino background 



density Qdm must reside in one or more species with 
a negligible (thermal) velocity to allow for structure 
formation. Moreover, relying on the astrophysical and 
cosmological considerations that point to the existence 
of dark matter (see p|fTO] and references therein), we 
think that a particle physics candidate for dark matter 
has to be electrically neutral, color neutral, and stable 
or have a lifetime Tdm that is generally larger than the 
age of the Universe today, to ~ 14 Gyr, i.e.Q 



Tdm > to - 4.3 X 10" s 



(2) 



With the standard active neutrinos being too light, 
such a dark matter candidate cannot be found within 
the Standard Model. Thus, one can consider the exis- 
tence of dark matter as evidence for new physics. 
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^ Studies of diffuse x-ray and 7-ray backgrounds, of the cos- 
mic ionization history, and of the CMB point to radiative life- 
times of electromagnetically decaying dark matter that are or- 
ders of magnitude above to [lim2, 13, 14, 15|. Moreover, for 
dark matter with decays into weakly-interacting relativistic 
particles such as neutrinos, Tdm S> to is inferred from studies of 

of type la supernovae (SN) 



LSS, of the CMB 
and of contraints on the cosmic neutrino fiux 
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In this review we focus on dark matter candidates 
that appear once the Standard Model is extended with 
the Peccei-Quinn (PQ) symmetry and/or supersymmetry 
(SUSY): the axion, the lightest neutralino, the gravitino, 
and the axino. These hypothetical particles are particu- 
larly well motivated: 

— The PQ mechanism allows for an elegant solution of 
the strong CP problem [2Tll22| . In fact, an additional 
global U(l) symmetry — referred to as PQ symmetry — 
that is broken spontaneously at the PQ scale p3 l [24 l [2] 
/a ^ 6 X 10® GeV can explain the smallness (or van- 
ishing) of the CP violating 0-vacuum term in quan- 
tum chromodynamics (QCD). The associated pseudo- 
Nambu-Goldstone boson is the axion [251126] which be- 
comes massive due to QCD instanton effects. Indeed, 
the axion is electrically neutral, color neutral, and suf- 
ficiently long-lived for being a compelling dark matter 
candidate 27,28,29j. 

— SUSY extensions of the Standard Model are an appeal- 
ing concept because of their remarkable properties, for 
example, with respect to gauge coupling unification, 
the hierarchy problem, and the embedding of grav- 
ity [30 l l3T | l32 | l33 p 34 l [35] . As superpartners of the Stan- 
dard Model particles, new particles appear including 
fields that are electrically neutral and color neutral. 
Since they have not been detected at particle accel- 
erators, these sparticles must be heavy or extremely 
weakly interacting. Moreover, because of the non-ob- 
servation of reactions that violate lepton number L or 
baryon number i?, it is often assumed — as also in this 
review — that SUSY theories respect the multiplicative 
quantum number 

R = [-lfB+L+2S ^ (3) 

known as R-parity, with S denoting the spin. Since 
Standard Model particles and superpartners carry re- 
spectively even (+1) and odd (—1) R-parity, its con- 
servation implies that superpartners can only be pro- 
duced or annihilated in pairs and that the lightest su- 
persymmetric particle (LSP) cannot decay even if it is 
heavier than most (or all) of the Standard Model par- 
ticles0 An electrically neutral and color neutral LSP 
candidate — such as the lightest neutralino, the grav- 
itino, or the axino — can thus also be a compelling dark 
matter candidate. 

For each dark matter candidate X, it is crucial to cal- 
culate its relic density fix and to compare the result 
with the dark matter density i7dm, for which a nomi- 
nal 3(7 range can be inferred from measurements of the 
CMB anisotropics by the Wilkinson Microwave Aniso- 
tropy Probe (WMAP) satellite gT] 

nll^h^ = o.io5ig:°^S (4) 

with h = 0.73to;o3 denoting the Hubble constant in units 
of 100 kmMpc~^s^^. Note that the nominal 3(t range 

^ While R-parity conservation is assumed in this review, its 
violation is a realistic option; see, e.g., [36II37I[38II39II40| . 



is derived assuming a restrictive six-parameter "vanilla" 
model. A larger range is possible — even with additional 
data from other cosmological probes — if the fit is per- 
formed in the context of a more general model that in- 
cludes other physically motivated parameters such as non- 
zero neutrino masses [42]: 0.094 < Qdmh"^ < 0.136. More- 
over, there are limits on the present free- streaming veloc- 
ity Vq of the dark matter candidate X from observations 
and simulations of LSS. Indeed, as mentioned above, the 
dominant part of i7dm has to have a negligible (thermal) 
velocity to allow for structure formation. 

Primordial nucleosynthesis — or big-bang nucleosynthe- 
sis (BBN) — is another cosmological probe for the viabil- 
ity of dark matter scenarios and, more generally, for the 
viability of models beyond the Standard Model. Rely- 
ing on Standard Model physics, general relativity, and 
the baryon-to-photon ratio inferred from the CMB an- 
isotropics, standard BBN (SBBN) predicts the primordial 
abundances of deuterium, helium, and lithium in good 
overall agreement with the observationally inferred val- 
ues and provides thereby an important consistency check 
of standard cosmology [UlliSliillTliHli] . This agree- 
ment can be affected by new physics, e.g., in the three 
following ways. First, there can be changes in the tim- 
ing of the nuclear reactions — caused, e.g., by a change in 
the Hubble expansion rate because of an increase of the 
energy density during BBN contributed by a new relativis- 
tic species [Mli451i50] . Second, non-thermal processes — 
caused, e.g., by the injection of energetic Standard Model 
particles in late decays of heavier particles into an ex- 
tremely weakly interacting dark matter candidate — can 
reprocess the produced light nuclei 51, 52 .53 , 54. 55.,56 ] [57 ] 
[58 l l59 l [60 l l6 1 (1621163] . Third, long-lived electromagnetically 
or strongly interacting relics — which can occur, e.g., in sce- 
narios with an extremely weakly interacting dark matter 
candidate — can form bound states with light nuclei and 
thereby lead to a catalysis of nuclear reactions, i.e., to 

\7 WWW7UWi^Wr9i . We will encounter the associated 
BBN constraints and other cosmological and astrophys- 
ical constraints explicitly in the discussions of the dark 
matter scenarios given below. 

The viability of a dark matter scenario is tightly con- 
nected to the early history of the Universe. As suggested 
by the flatness, isotropy, and homogeneity of the Universe, 
we assume that its earliest moments were governed by in- 
flation J80','81] . The inflationary expansion is then followed 
by a phase in which the Universe is reheated. The reheat- 
ing phase repopulates the Universe and provides the initial 
conditions for the subsequent radiation-dominated epoch. 
We refer to the reheating temperature Tr as the initial 
temperature of this early radiation-dominated epoch of 
our Universe. In fact, inflation models can point to 
well above 10^° GeV [HCTIICTS^ . Here one should stress 
that BBN requires a minimum temperature of [83U84ll85] 
[86] 

T> 0.7-4 MeV EE r„iin (5) 

and that any temperature above Tmin is still speculative 
since BBN is currently the deepest reliable probe of the 
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early Universe. Interestingly, futuristic space-based gra- 
vitational-wave detectors such as the Big Bang Observer 
(BBO) or the Deci-hertz Interferometer Gravitational Wa- 
ve Observatory (DECIGO) [87^ could allow for tests of 
the thermal history before BBN and could even probe the 
reheating temperature Tr after inflation [881189] . 

Let us now turn to the dark matter candidates one 
by one: axions (Sect. [2]), neutralinos (Sect. [3]), graviti- 
nos (Sect. |4|), and axinos (Sect. [5]). For each of those 
candidates, I review primordial production mechanisms, 
cosmological/astrophysical constraints, and prospects of 
experimental identification. 



2 Axion Dark Matter 

In this section we consider the axion which is a promising 
dark matter candidate that is not tied to SUSY but that 
can coexist with any SUSY dark matter candidate. Let 
us start this section with a brief review of the strong CP 
problem, the PQ mechanism, and a short description of 
basic properties of the axion. More detailed discussions 
can be found, e.g., in [8T 1 l90 l [9T l l24 l l2 l l92] . 

The isoscalar r( meson being too heavy to qualify as a 
pseudo-Nambu-Goldstone boson of a spontaneously bro- 
ken axial U(1)a is the well-known U(1)a problem of the 
strong interactions. This problem has an elegant solution 
in QCD based on non-trivial topological properties [931 
(94) . The r\' meson mass can be understood as a conse- 
quence of the Adler-Bell-Jackiw anomaly t95y96j receiving 
contributions from gauge-field configurations with non- 
zero topological charge such as instantons [97] . The ex- 
istence of such configurations implies the additional Q- 
vacuum term in the QCD Lagrangian [93l[94l[95l[96] 



(6) 



with the strong coupling and the gluon-field-strength 
tensor GJ^j, whose dual is given by G^^ = e^^paG°''"^/2. 
The ©-vacuum term violates the discrete symmetries P, 
T, and CP for any value oi ^ mi with n e Z. Such vio- 
lations have not been observed in strong interactions and 
experiments on the electric dipole moment of the neutron 
give an upper bound of \0\ < lO'^. Within QCD, = 
seems natural based on the observed conservation of those 
discrete symmetries. However, once QCD is embedded in 
the Standard Model of strong and electroweak interac- 
tions with CP violation being an experimental reality, 
and the argument of the determinant of the quark mass 
matrix M — two a priori unrelated quantities — must can- 
cel to an accuracy of 10"^ according to the upper bound[3 



\0\ = |6> + ArgdetMl < 10" 



(7) 



This fine-tuning problem is the strong CP problem (cf. 
also [90ll9T] and references therein). 



^ If one or more quarks are massless 
away by a chiral rotation; cf. [901191] . 



can be rotated 



The elegant solution of the strong CP problem sug- 
gested by Peccei and Quinn [^[^ requires a new global 
chiral U(l) symmetry — the PQ symmetry U(1)pq — that 
is broken spontaneously at the PQ scale fa- The cor- 
responding pseudo-Nambu-Goldstone boson is the axion 
a [25ll26] . which couples to gluons such that the chiral 
anomaly in the U(1)pq current is reproduced. 



~ fa/N 32^2 ^M-^ 



(8) 



where the color anomaly of the PQ symmetry depends 
on the axion model as discussed below. This interaction 
term ^ together with the vacuum term ([6]) for — > 
= 0-1- Arg det M provide the axion field with an effec- 
tive potential VcS at low energies. This solves the strong 
CP problem since the coefficient of the CP violating GG 
term becomes dynamical and vanishes for the value (a) = 
~0/a/A^ at which VcH has its minimum. 

Because of the chiral U(1)pq anomaly, the axion re- 
ceives a mass from QCD instanton effects [MlllOO| that 
govern the effective axion potential V^s at low energies 



ruumd 



f-^rn-j^ 
(m„ + md)2 \fa/N 



(9) 



where (mj) is the mass of the up (down) quark and 
« 92 MeV and = 135 MeV are respectively the 
decay constant and mass of the pion. Using z = rriu/md, 



ma 



1 + 2 \fa/N 



0.60 meV 



2 /lOi°GeV\ 
(iTo^V fa/N ) 



1 



(10) 
(11) 



This relation is shown for z — 0.56 [101] and fa/N ^ fain 
Fig. □ (from [25); note that z could be within 0.3-0.6 [1]. 

While the original PQ proposal [^Tll^ assumed fa to 
be at the weak scale, axion searches, astrophysical obser- 
vations, and cosmological arguments, which are discussed 
below, point to a significantly higher value of the PQ scale 
(cf. [23u24u 2j and references therein) 



fa/N > 6 X 10^ GeV 



(12) 



Accordingly, the interactions of the axion are strongly sup- 
pressed and its mass must be very small, rria 0.01 eV, 
so that the axion can be classified as an extremely weakly 
interacting particle (EWIP). 

Axion interactions are model dependent. The two most 
popular classes of phenomenologically viable "invisible ax- 
ion" models are the hadronic or Kim-Shifman-Vainshtein- 
Zakharov (KSVZ) models [l02l[T03] and the Dine-Fischler- 
Srednicki-Zhitnitskn (DFSZ) models [I04llT05] . In models 
of the KSVZ type, at least one additional heavy quark is 
introduced which couples directly to the axion while all 
other fields do not carry PQ charge. Thus, the axion in- 
teracts with ordinary matter through the anomaly term 
from loops of this new heavy quark. Integrating out the 
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Fig. 1. The PQ scale fa, its relation to the axion mass ma 
(for A'^ = 1) together with axion-search ranges (dark-shaded) 
and exclusion limits from axion searches, cosmological con- 
straints, and astrophysical limits (light-shaded and medium- 
shaded, columns from left to right). The vertical line with the 
two arrowheads indicates the fa range where Oa from the mis- 
alignment mechanism can provide naturally a sizable fraction 
of J7dm for fa > Tr. For a scenario with fa < Tr, Oa > ^?dm in 
the region shown by the "Cold DM" bar. The other bar in this 
column indicates the region that is excluded by constraints 
from LSS ("Hot DM") or by a thermal relic axion density 
^therm ^ ^^^^ The trausitiou region between those constraints 
is excluded by telescope searches for a — > 77 decays as indi- 
cated by the "Telescope" bar in the leftmost column. The dark- 
shaded bars in that column indicate the search ranges of the 
axion-dark-matter-search experiment ADMX and of the axion 
helioscope CAST. The exclusion ranges in the rightmost col- 
umn are inferred from the supernova SN 1987A, i.e., from the 
observed duration of the emitted neutrino burst ( "Burst dura- 
tion") and from the non-observation of axion-induced events 
("Too many events"). The second column from the right ap- 
plies to models with a direct axion-electron coupling and is 
inferred from globular cluster (GC) stars and white-dwarf cool- 
ing. Also the middle column is inferred from GC stars but from 
the axion-photon coupling: While the medium-shaded bar ap- 
plies to a DFSZ model, the exclusion range can be more ex- 
tended in a KSVZ model as indicated by the light-shaded tip. 
Details are given in the main text. From [24j . 



heavy quark loops, one obtains the effective dimension- 
5 couphng of axions to gluons given in ((8|) with A'^ = 1. 
In particular, couplings of the axion to Standard Model 
matter fields are suppressed by additional loop factors. 
In the DFSZ schemes, no additional heavy quarks are in- 
troduced. Instead, the Standard Model matter fields and 
at least two Higgs doublets carry appropriate PQ charges 
such that the axion also couples directly to the Standard 
Model fields. Again, at low energies the axion-gluon inter- 
action ([8]) arises, but now with N = 6. For more details 
on these axion models, we refer to the reviews [106, 107[ 
[92] . In addition, axions with fa ~ 10^^ GeV appear also 
in string theory; cf. p ^ ll09 p il0j. 

For the axion lifetime Ta and axion phenomenology, 
the axion-photon interaction plays a central role: 



r 



9a 



77 p T^Hf 



(13) 



with the electromagnetic field-strength tensor F^^ , its dual 
F^i, = and the coupling constant 



E 24 + z 



a 1 + Z THa 

27r y/z f^m^ 



31 + z 
E 24 + z 
N ^ 31+z 



(14) 
(15) 



where both the electromagnetic anomaly E and the color 
anomaly N of the PQ symmetry depend on the axion 
model. For example, E/N = in a KSVZ model with 
(an) electrically neutral new heavy quark(s) and E/N = 
8/3 in DSFZ models. Note that E/N = 2 is a possibil- 
ity for which the axion-photon coupling would be sup- 
pressed [111) due to cancellations in the bracket of (114115^ . 

The lifetime of the axion governed by the decay a — > 77 
(and the associated rate Fa^^j) is obtained from P^ : 



a— >77 



4.6 X 10*° s 



647r 

E 
N 



- 1.95 



fa/N Y 
1010 GeV / 



(16) 



(17) 



where z = 0.56 was used. Comparing Tq for E/N = 
with the present age of the Universe to given in ^ , one 
finds Ta > to for fa/N > 3 X 10^ GeV or equivalently 
nia ^ 20 eV. Thus, the axion is stable on cosmological 
time scales for the fa/N values that are favored by 
cosmological and astrophysical constraints. 



2.1 Primordial Origin 

Axion dark matter can have different origins and different 
properties depending on the cosmic history and on the 
value of fa/N. A remarkable generic property of the axion 
is the fact that it is massless for T > 1 GeV > ^qcd and 
that it acquires mass only through instanton effects for 
T < Aqcd such that ma is given by for T ^ 0. 

The most straightforward situation is encountered when 
fa/N is sufficiently small so that the axion couples suffi- 
ciently strong for being in thermal equilibrium with the 
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early primordial plasma. The axion will then decouple (or 
freeze out) at a temperature Tf 3> rua as a thermal relic 

0^077 (^)(^) (18) 

where g^,s denotes the number of effectively massless de- 
grees of freedom such that the entropy density reads s — 
2TT'^g^sT^/^5. For fa/N < 3 x 10^ GeV (corresponding to 
rUa ^ 0.2 eV), the axion is a thermal relic that decou- 
ples after the quark-hadron transition, Tf < 150 MeV, 
where tttt <-!■ tto is a generic process that keeps the axions 
in thermal equilibrium before their decoupling [112illl3j 
Accordingly, the dark matter constraint f2a ^dm implies 
iTia ^ 18 eV, where also > to is satisfied as discussed 
above. However, such axions are hot dark matter so that 
f^a can only be a minor fraction of J^dm- In fact, obser- 
vations of LSS imply a restrictive limit of ma ^ 1-02 eV 
(95% CL) 113,7,116J, which is indicated as "Hot DM" in 
Fig. [1] (second column from the left). 

For large values of fa/N such that axions are never in 
thermal equilibrium with the primordial plasma, f2a be- 
comes sensitive to the earliest cosmological history. In this 
limit, the simplest setting is encountered when the sponta- 
neous breaking of the PQ symmetry occurs before inflation 
leading to a reheating temperature of Tr < fa so that no 
PQ symmetry restoration takes place during inflation or 
in the course of reheating. In this setting, axionic strings 
can be neglected since they have been diluted by inflation 
and cannot be produced later on. Moreover, the initial 
misalignment angle &i of the axion with respect to the 
CP-conserving position (i.e., its position in the mexican 
hat potential) is everywhere (basically) the same in our 
observable patch of the Universe. Since Oi is typically not 
at the minimum of the effective potential (i.e., not at the 
minimum of the tilted mexican hat) generated by instan- 
ton effects at T ~ ^qcd , the field starts to oscillate coher- 
ently around its minimum for ma(Tosc) — 3H{Tosc)- Once 
rUa takes on its T-independent value (jQllOp . this axion 
condensate behaves as cold dark matter 27,28,29 with a 
relic density that is governed by the initial misalignment 
angle -n < O, < n [TTTIITTS] : 

^ah^^O.15U{0l)O!l^j^^y' (19) 

with ^ = 0(1) parametrizing theoretical uncertainties re- 
lated, e.g., to details of the quark-hadron transition and of 
the T-dependence of rUa (cf. |118j ) and f{Of) accounting 
for anharmonicity at sizable 6*^; f{0^) — > 1 for 6>f 0. 
For lO^'^GeV < fa/N < lO^^ Qgy, this "misalignment 
mechanism" can provide naturally, i.e., for 0i = 0(1), a 
sizable contribution f2a to f^dm, which is indicated by the 
vertical line with the two arrowheads in Fig. [1] However, 
with any value of 16*^1 G [0,7r] being equally probable, 
there is always the possibility of 16*^1 « in our patch 

* Before the quark-hadron transition, axions can be kept in 
thermal equilibrium by QCD reactions that involve the axion- 
gluon interaction ^ with a sufficiently small fa/N [11411115] . 



of the UniverseH In fact, for a finely tuned |63,;| < 10"^, 
this setting allows for the high values of fa ~ 10^^ GeV 
suggested by concepts of grand unification and string the- 
ory |108pi09pil0j . With the original motivation for the 
axion residing in a solution of the fine-tuning problem ^ , 
a finely tuned 0i is somewhat unpleasant but can be asso- 
ciated with anthropic considerations |119pi20 j. Note that 
it would actually be more correct to replace Of in by 
(6*^) at r ~ ^QCD, which cannot be arbitrarily small in 
the case of infiation with energy scale Ei or Hubble scale 
Hi — SnG-^Ef /3 defined when those modes exit the hori- 
zon that reenter the horizon today p ^ ll22 p i23 p i24 p i08 | 

ma, 

due to fluctuations in the axion field from de Sitter quan- 
tum fiuctuations during the inflationary epoch. 

Let us now consider the case in which the spontaneous 
breaking of the PQ symmetry occurs after infiation, Tr > 
fa, again for large values of fa/N such that axions are 
never in thermal equilibrium with the primordial plasma. 
After spontaneous breaking of the PQ symmetry at T ~ 
fa, our observable patch of the Universe will consist of 
many smaller patches, each of which with an arbitrary 
value of Oi G (— tt, tt], and of an associated network of 
axionic strings and domain walls. Indeed, since a uniform 
distribution of Oi is provided with (Of) = 7r^/3, the axion 
relic density from the misalignment mechanism becomes 
independent of Oi : 

which includes the anharmonic correction factor /(Q™^*^ = 
tt/VS) — 1.2. The dark matter constraint Ha < i7dm thus 
implies, e.g., for a moderate ^ ~ 0.4 the limit fa/N < 
6x 10-^^ GeV (ma > 10~^ eV) as indicated in Fig. [1] (second 
column from the left). This limit is conservative since there 
can be additional sizable contributions to f2a from decays 
of axionic strings, domain walls, and non-zero momentum 
modes of the axion field. For a discussion of those contri- 
butions and potential domain wall problems, see [115] and 
references therein. 

2.2 Cosmological Constraints 

The thermal relic axions encountered for sufficiently small 
values of fa/N will contribute to the radiation density at 
the time of BBN in a way analogous to the case of an extra 
light neutrino species. Thereby, these axions will give the 
following contribution to the effective number of neutrinos 
Nf at T ~ 1 MeV dH] 

^ In fact, as a matter of principle, it is impossible to calculate 
the precise value of Sla from first principles in this setting. 
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For Tt > 1 MeV, ANf < 4/7 = 0.57 which is compatible 
with the BBN hmit |l9| 

Nf = 3.lt\i (95% CL) ; (23) 

see also Fig. 2 of Ref. |;50^. However, in the literature also 
more restrictive limits on AN^^ can be found, which can 
be associated with an nia bound that is more restrictive 
than the limit nia < 1.02 eV (95% CL) pl3 l [7 l lTT6 ] from 
observations of cosmological LSS discussed in the previous 
section. For example, the limit AN^^ < 0.25 would imply 
g^s{Ti) > 20 and thereby fa > 2 x 10^ GeV (or < 
0.3 eV) as can be seen from Fig. 4 of Ref. [113j . 

For large fa/N scenarios with PQ breaking before in- 
flation and Tr < fa, in which axions are never in ther- 
mal equilibrium with the primordial plasma, axion density 
fluctuations produced during inflation lead to isocurvature 
[125,126,127,122,82 ,117,120] and non-Gaussian 123,12.8, 
1129] effects in the temperature fluctuations of the CMB. 
Such effects are constrained by cosmological data. 

With the axion field being essentially massless during 
inflation, de Sitter quantum fluctuations are imprinted on 
the axion as on any other light scalar field. Since the cou- 
plings of the axion are so weak, an isocurvature mode 
survives that is uncorrelated to the usual adiabatic mode 
seeded by the quantum fluctuations of the inflaton field. 
The resulting deviations from adiabaticity could manifest 
themselves in the temperature fluctuations of the CMB 
and could thereby support the existence of the axion (or, 
more generally, of at least two light scalar fields during 
inflation, one of which giving rise to dark matter that is 
never in thermal equilibrium with the primordial plasma) . 

So far, no evidence for dark matter isocurvature per- 
turbations has been found. Using ao to parameterize the 
ratio of the power spectrum of entropy perturbation, Ps(fc), 
to the one of curvature perturbation, Pji{k), at the pivot 
wavenumber fco — 0.002 Mpc~^, 

ajko) ^ Ps jko) ^24^ 
1 - a{ko) - Pniko) ' ^ ' 

current limits on axionic entropy perturbations are ao < 
0.16 (95% CL) and ao < 0.067 (95% CL) obtained re- 
spectively from a WMAP-only and a WMAP-hBAO-l-S]N0 
analysis [Tj. Confronting the ao limit with the expression 
derived for an axion that receives quantum fluctuations 
with a nearly scale invariant spectrum during inflation 
and TV = 1 [mjH 

„,„...4xl0-2^, (25) 

one obtains cosmological constraints on fa that depend on 
El (or Hi) and on the fraction Ra = /?a/^?dm [TTTlfnUj . 
Figure[2](from |120j) shows isocurvature limits inferred 

® Here BAO refers to baryon acoustic oscillations and SN to 
supernovae of type la. See [l] for references and details. 

Note that there are typos in Eqs. (44) and (27) of Ref. PT7] . 
The corrected prefactor in (44) is 2.4 x 10^° and the associated 
limit in (27) should read 1.2 x 10^** GeV. I thank J. Lesgourgues 
for the clarification of this point. 



Hubble Lcalc of inflation H; (GcV) 
10« IQP 10'° 10" 10" 10'^ 10'" 10'= 




10'^ 10'" 10'= 10"^ 10" 

Energy scale of Inflation £^7 (GeV) 



Fig. 2. Constraints on the PQ scale /„ (for A'^ = 1) as a 
function of the energy scale of inflation Ei or Hubble scale 
Hi = SnGr^Ef /3 defined when those modes exit the hori- 
zon that reenter the horizon today. The thick solid (red) lines 
show fa = Hi/2n (upper panel) and fa = Tmax = eoff-Ei for 
Eoff ~ 1, 10"'^'^, 10"'' (lower panel). Below (above) the cor- 
responding line, PQ symmetry restoration does (not) occur 
during/after inflation. The hatched regions above those lines 
show the cosmological constraints from the WMAP-f-BAO-|-SN 
limit ao < 0.067 (95% CL) on axionic isocurvature pertur- 
bations (|25|) . Those constraints depend on Ra = i7a/f?dm as 
labeled and as indicated by the different shadings. The con- 
straints labeled as "Any Ra" result from the lower limit on 
axion fluctuations due to quantum fluctuations during infla- 
tion ((20I). For El > 10^'^ GeV, this can lead to Oa > J7dm as in- 
dicated by the (green) region labeled as "Too much CDM" . The 
WMAP-(-BAO-|-SN limit on the tensor to scalar ratio r < 0.2 
(95% CL) excludes the region with Ei > 3.8 x 10^^ GeV. 
The exclusion ranges from axion searches and astrophysical 
considerations are indicated respectively by the dark (brown) 
horizontal bar and by the light-shaded (yellow) region with 
fa < 10^ GeV. For fa > Mp, the PQ scale is super-Planckian. 
More details are given in the main text. From '120'. 
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from ao < 0.067 for Ra — I (cyan) and Ra > 0.25 (purple) 
and f2a < ^?dm limits (green). The Ra independent limits 
(pink) are associated with the lower limit ([20]) . For each of 
the constraints that depend on p9|) . a more and a less con- 
servative limit obtained respectively for ^ = 0.05 and ^ = 
1 is indicated by the darker and lighter shadings. The kink 
of these limits for fa > 10^^ GeV results from Tqsc ^ ^qcd 
and an associated Ha expression with a different paramet- 
ric dependence on fa/N than in (fTQ]) derived for Tqsc 
^qcd; for details, see e.g. [108pi20j . The upper (lower) 
panel shows the case of very inefficient (efficient) thermal- 
ization at the end of inflation. The thick solid (red) lines 
indicate fa = fl"i/27r and fa = T^ax = Coff^^i for eoff = 1, 
10~^'^, 10"'^. Below (above) the corresponding line, PQ 
symmetry restoration does (not) occur during/after in- 
flation. Obtained from a WMAP-I-BAO+SN analysis Q^, 
the current limit on the amplitude of primordial gravi- 
tational waves expressed in terms of the tensor to scalar 
ratio r < 0.2 (95% CL) excludes Ei > 3.8 x lO^^ GeV as 
indicated by the orange region. In the region indicated by 
the horizontal (brown) bar, axion models with an unsup- 
pressed axion-photon coupling ga-y-y are excluded by the 
axion dark matter search experiment ADMX; cf. Fig. [3] 
and discussions in the next section. The (yellow) region 
with fa < 10^ GeV is disfavored by conservative astro- 
physical constraints; cf. Fig. [T] The dashed lines indicate 
targets of future experiments/observations. 

While the quantum fluctuations in 0i during infla- 
tion are Gaussian distributed, the axion-induced temper- 
ature fluctuations in the CMB are governed by fluctua- 
tions in the axion number density after formation of the 
axion condensate with Ua (x Of modulo anharmonic cor- 
rections. Thereby a non-linearity enters that leads to a 
non-Gaussian contribution to the CMB temperature fluc- 
tuation 123, 128, 129). This may allow one to probe axions 
cosmologically even for the case fla ^ ■!^dm [129] . 

If the PQ symmetry is not broken before inflation, 
there will be no light axion field present to acquire de 
Sitter quantum fluctuations and thus no axionic source 
of isocurvature perturbations. Also if the PQ symmetry 
is restored during inflation, Hi/2-k > fa (i.e., the typical 
amplitude of quantum fluctuations exceeds the PQ scale), 
or during reheating, Tmax > fa (i-e., thermal PQ sym- 
metry restoration), primordial axionic perturbations will 
be washed out which implies negligible axion-induced de- 
viations from adiabaticity. Moreover, as described in the 
previous section, the subsequent PQ breaking will lead to 
a sample of Oi values with a flat probability distribution 
after PQ breaking. Since these Oi values are distributed at 
random throughout our observable Universe, spatial axion 
density variations are expected once the axion acquires its 
mass. These variations can then lead to "axion mini clus- 
ters" containing a sizable contribution to J7dm [115] . 

2.3 Astrophysical Constraints 

Axions could be produced not only in the early Universe 
but also in hot and dense astrophysical environments such 
as those encountered in ordinary stars, white dwarfs, and 



supernovae. The axion luminosity La of such sources de- 
pends on /a, the relevant axion production processes (and 
thereby on the axion couplings and model) , and on the as- 
trophysical understanding/model of the source under con- 
sideration. A sizable axion luminosity is associated with an 
additional energy transport out of the corresponding as- 
trophysical source. This can affect the behavior/character- 
istics of the source strongly. Astrophysical studies of stars, 
white dwarfs, and supernovae can thus be used to derive 
constraints on fa/N or ma [l30lll31lll32l[T33l[8Tl[T34l[24] . 

To allow for axion production, the relevant energy/tem- 
perature in the source must be sufficiently large with re- 
spect to TOa already for kinematical reasons. In addition, 
the interaction strength of the axion enters in the following 
two ways: On the one hand, a stronger axion interaction 
(smaller fa/N or larger rua) allows for a more efficient pro- 
duction in the source. On the other hand, stronger inter- 
actions are also associated with smaller mean free paths in 
medium and possibly with rescatterings within the source 
which can delay energy emission via the axion channel. 

The current status of astrophysical axion bounds is 
reviewed, e.g., in Refs. p4l[2] and shown in Fig. [1] (three 
rightmost columns). Here only a brief summary is given. 

In stars including our sun, axions can be produced 
via their coupling to photons ga-y-y through the Primakoff 
process j + Ze -> Ze + a [l30 l ll35 p i36 ( fT37]. The ax- 
ionic energy drain (described by La (x (7^77) can then 
lead to an enhanced consumption of nuclear fuel within 
the star and can thereby, e.g., shorten the lifetime of a 
star. In this respect globular clusters (GCs), which are 
bound systems of a homogeneous population of low-mass 
stars, are particularly valuable since they allow for tests 
of stellar-evolution theory. In fact, studies of GC stars 
point to 5a77 < 10"^° GeV"^ 138,134,24J which implies 
fa > 2.3 X 10^ GeV in a KSVZ axion model with E/N = 
and fa > 0.8 x 10'' GeV in a DFSZ axion model with 
E/N — 8/3 for z = 0.56. These limits are indicated by 
the middle column in Fig. [TJ Not shown are the limits 
inferred from studies of our sun since they are less restric- 
tive. For example, helioseismological studies of the sound- 
speed profile give the limit ga-yy < 10~® GeV~^ |139j . 
Moreover, ga^-y ^ 5 x 10"^" GeV"'^ [M] is inferred from 
measurements of the solar ^B neutrino flux which is sen- 
sitive to enhanced nuclear burning [139j . 

Axion production via 7 -f e~ — > e~ -I- a and e~ + Ze 
Ze+e~ +a is more efiicient than the Primakoff processes in 
models with a direct axion-electron coupling gaee- Studies 
of GC stars [140) and of white-dwarf cooling |14HI142[ 
I143[|144j thereby lead to restrictive limits on gaee- For the 
DSFZ model with two Higgs doublets carrying PQ charge 
and the ratio of the associated Higgs vacuum expectation 
values (VEVs) given by tan/3, these fimits imply fa < 
1.3 x 10^ GeV cos^ /3 [24J which is indicated for cos^ /3 = 0.5 
by the second column from the right in Fig. [TJ 

In core collapse supernovae (SN of type II) that lead to 
the formation of a hot — T = 0(10 MeV) — proto-neutron 
star, axions can be produced via their coupling to nucleons 
QaNN through axion bremsstrahlung emission N -\- N ^ 
N + N + a [145pi46) in the dense nuclear medium [147j 



8 



Frank Daniel Steffen: Dark Matter Candidates 



I145pi48pi49pi50pi51) . In fact, in these dense environments 
even the mean free path of neutrinos is such that they 
rescatter and diffuse out before carrying away energy. This 
picture was confirmed by the observed duration of about 
10 s of the neutrino burst from the supernova SN 1987A. 
Since axion emission would be an additional energy drain, 
it would reduce the cooling time and this burst dura- 
tion. The observed burst duration thereby implies fa ^ 
4 X 10* GeV [24] as indicated by the upper bar in the 
rightmost column of Fig. [T] For /a < 6 x 10^ GeV, ax- 
ions cannot affect the neutrino burst duration significantly 
since their coupling becomes so strong that they rescatter 
within the dense medium [T521IT53]. For /a < 3 x 10^ GeV, 
the axion coupling would have been strong enough for ax- 
ion induced events that have not been observed for the 
supernova SN 1987A 1151] . This is indicated by the lower 
bar in the rightmost column of Fig.[TJ The gap between the 
two bars is known as the "hadronic axion window" which 
was thought to be an allowed region for KSVZ models 
with strongly suppressed ga^y-y [155] . However, this win- 
dow is closed by the hot dark matter constraint from LSS 
observations |113P 71 [TTB] discussed in Sect. 12.11 and indi- 
cated by the second column from the left in Fig. [TJ 



2.4 Experimental Searches and Prospects 

Axion searches started already more than 30 years ago and 
excluded soon fa values close to the weak scale, e.g., from 
studies of the branching ratio BK{K^ — + 7r+ -I- nothing), 
as indicated by the bar at the bottom ("Laboratory") of 
the leftmost column in Fig. [TJ Present day axion searches 
are probing much larger values of fa/N including those 
of invisible axion models and those in which axions can 
provide the dominant component of dark matter ^23iil56i, 
[2]. Depending on the expected origin of the axions, these 
searches can be classified into those for cosmic axions, so- 
lar axions (or more generally axions from astrophysical 
sources), and axions produced in the laboratory. 

Most axion searches rely on the axion-photon coupling 
5077 Since the axion has not been discovered so far, neg- 
ative searches can be translated into ma-dependent (or /q- 
dependent) limits on ga-yy ■ A summary of those exclusion 
limits is given in Fig. [3] which is an update of Fig. 10.26 
of Ref. [156j by courtesy of M. Kuster. Here the diagonal 
(yellow) bar and the solid diagonal line indicate the region 
of typical invisible axion models and gaj-y in the KSVZ 
model with E/N — 0, respectively. The cosmological hot 
dark matter constraint nia ^ 1 eV |113[ [7l fTT6] is shown 
by the short vertical (orange) line ("HDM") within this 
region and the astrophysical limit from studies of stars on 
the horizontal branch (HB) in GBs gajj < 10"^° GeV"^ 
by the vertical dashed (orange) line ("HB stars") [138[ 
11341 124] . Let us recall that axions can provide a significant 
contribution to f^dm for nia ^ 3 x 10~^ eV. 

Sizable cosmic axion densities are expected in galaxy 
clusters and galaxies and thereby also on Earth. There are 

® Fifth-force experiments allow for ga-^^^-independent axion 
searches |157l : see [158] and references therein. 




lO"'" 10"" 10"' lO"-* 10"' 10° 10^ 



Fig. 3. Axion exclusion limits in the plane spanned by ma 
and ga-yy - The diagonal (yellow) bar ("Axion Models") and the 
solid diagonal line ("KSVZ E/N=0") show respectively the re- 
gion of typical invisible axion models and gaj-y in the KSVZ 
model with E/N = 0. The cosmological hot dark matter con- 
straint [113| [7l fTT6| is shown by the short vertical (orange) line 
("HDM") within this region and the astrophysical limit from 
studies of stars on the horizontal branch (HB) in GBs by the 
vertical dashed (orange) line ("HB stars") p8lll34U24] . The 
telescope search for a — > 77 decays reported in Ref. [159) ex- 
cludes the pencil-like (orange) region ( "Telescope" ) and direct 
axion searches with microwave cavities the pencil-like gray re- 
gion ( "Microwave Cavity" ) [160lll6llll56[ . Solar axion searches 
with Eixion helioscopes exclude the regions above the upper 
("Lazarus et al") and lower ("Tokyo Helioscope") dotted lines 
and the one above the solid (blue) line ("CAST") [T62lll63U164l 
I165| . The medium-shaded (cyan) region ("Bragg Reflection") 
is excluded by searches for axion conversion in the Coulomb 
field of nuclei in a crystal lattice [1661116711168] and the light- 
shaded (green) region ( "Laser Experiments" ) by searches for a 
"light-shining through a wall" event [169] . I thank M. Kuster 
for providing this update of Fig. 10.26 of Ref. [156] . 



indirect and direct searches for those cosmic axions. In the 
indirect ones, telescopes are used to look for photons from 
a ^ 77 decays, for example, in the Abell clusters |170pi711 
1159) and in nearby dwarf galaxies [172] . After correcting 
for the Doppler shift due to the motion of the host, one 
expects a basically mono-chromatic spectrum of the re- 
sulting 7's at — ma/2. In addition to possible insights 
into nia, also is probed in these telescope searches. No 
such 7's have been identified unambiguously so far which 
implies that Ta can respect ([2]). The associated exclusion 
limits (labeled as "Telescope") are shown by the corre- 
sponding bar in the leftmost column in Fig. [1] and by the 
vertical pencil-like (orange) region in Fig. [S] 
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In direct searches for cosmic axions, microwave cav- 
ities are used to fook for resonant conversion of those 
axions — pervading the Earth — into photons |173ill74ill75[ 
I176pi60|ll6ipi77j afong the hafoscope technique proposed 
in |178j . Through the cavity frequency at which the res- 
onance would appear and throught the width of the res- 
onance, these experiments are sensitive to ma and to the 
virial distribution of thermalized axions and thereby to 
the axion distribution in the galactic halo. In fact, mi- 
crowave resonant cavity experiments probe exactly the 
nia range in which axions can provide a sizable contribu- 
tion to S^dm', cf. bar labeled as "ADMX" in the leftmost 
column in Fig. [TJ Moreover, the ADMX experiment has 
achieved a sensitivity such that realistic axion models have 
already been probed and excluded at 90% CL in a narrow 
nia range [16011161) . Indeed, no axion signal has been ob- 
served so far. The associated exclusion limits are shown 
by the horizontal (brown) bar in Fig.[2]and by the vertical 
pencil-like gray regions labeled as "Microwave Cavity" in 
Fig. [3l An upgrade of ADMX is underway that should al- 
low one to probe realistic axion models over a much larger 
nia range |1761ll77j . Relying on Rydberg-atom detectors, 
an upgrade of the cavity experiment CARRACK is aiming 
at a similar sensitivity and search range [179^ 176 . 

As already addressed in the previous section, axions 
could be produced in astrophysical sources such as our 
sun. Searching for solar axions means to look for the con- 
version of such axions into 7's in an electromagnetic field 
via the inverse Primakoff process. In axion helioscopes [1781 
fT80lfT8n[T82l[T62l[T63] . this field is provided by a strong 
magnet that is pointed at the sun. Using Bragg diffrac- 
tion at crystal detectors [168, 166, 167,, it is the Coulomb 
field of the nuclei in the crystal lattice that provides the 
field. Also the geomagnetic field of the Earth can allow in 
principle for the conversion of solar axions into photons 
which could be detected by a satellite on the dark side 
of the Earth |183| . Existing solar axion searches provide 
so far only rria-dependent exclusion limits on gaj-y- The 
ones from axion helioscopes are indicated by the upper 
("Lazarus et al.") and lower ("Tokyo Helioscope") dotted 
lines and by the sohd (blue) hne ("CAST") line in Fig. [3 
Among those limits, the most restrictive one is provided 
by the CERN Axion Solar Telescope (CAST): ga^^ < 
8.8 X 10-" GeV"^ for < 0.02 eV .162,, 163 . Currently, 
the Tokyo Helioscope and CAST are probing already the 
parameter region with realistic axion models |1641I165| but 
for nia > 10"^ eV which is associated with warm/hot 
axion dark matter scenarios. Indeed, an axion discovery 
in that region would imply that axions cannot be the 
dominant component of f2dm for a standard cosmologi- 
cal historjd so that there is room, for example, for the 
LSP to take over that part. The TOo-independent exclu- 
sion limit from Bragg reflection ga-yj ^ 1.7 x 10~^ GeV~^ 
obtained by the DAMA collaboration |167j is indicated by 
the medium-shaded (cyan) region (labeled accordingly) in 
Fig. [3l similar limits were obtained by COSME [166J and 
SOLAX [168]. Not shown is the limit gay^ < 10"" GeV 

^ For axion cosmology with a non-standard cosmological his- 
tory, see e.g. [184] . 



for nia S 10~^ eV that has been inferred from the absence 
of 7-ray bursts in coincidence with the SN 1987A neu- 
trino burst {1851 11861 . Such bursts could have originated 
from axions produced in the SN 1987A that had been con- 
verted subsequently into photons in the galactic magnetic 
field. 

Axion searches are also performed purely in the labora- 
tory. Using the Primakoff process, one should in principle 
be able to convert a small fraction of the photons of a laser 
beam in a strong transverse magnetic field into axions. In 
contrast to the photons, those axions should be able to 
traverse a wall. In a second strong magnetic field behind 
the wall, there should also be a non-zero probability for 
the inverse Primakoff effect in which the axion is recon- 
verted into a detectable photon. The overall probability of 
such a "light-shining through a wall" event is proportional 
to gaj-f- No event of this sort has been observed so far and 
the Brookhaven-Fermilab-Rutherford- Trieste (BERT) col- 
laboration has inferred ga-y-y < 6.7 x 10~'^GeV (95% CL) 
for nia < 10"^ eV 169 as indicated by the light-shaded 
(green) region labeled as "Laser Experiments" in Fig. [3] 

Another way to search for axions in the laboratory re- 
lies on the prediction that axions can affect the polariza- 
tion of light that propagates in vacuum through a trans- 
verse magnetic field. In fact, because of their coupling to 
photons (?a77, axions can induce dispersive and absorptive 
processes and thereby the following two phenomena [187j : 
(i) Linear dichroism, which refers to a rotation of the po- 
larization vector by a finite angle, and (ii) birefringence, 
which refers to the introduction of an ellipticity and an 
rotation of an initially linearly polarized beam. Searching 
for these phenomena, the BFRT experiment has extracted 
the exclusion limit 5^77 < 3.6 x lO"'^ GeV (95% CL) for 
nia < 5 X 10~^eV. Note that the evidence for vacuum 
dichroism claimed by the PVLAS collaboration in the 
year 2006 |188j has been retracted recently [189]. Origi- 
nally, they had interpreted their findings in terms of the 
presence of an axion-like particle (ALP) with a mass of 
(1—1.5) X 10~^ eV and a coupling to photons in the range 
(1.7-5) X 10^6 GeV"^ [HH]. The retraction reassures the 
validity of the astrophysical constraints that had already 
excluded the region in which this signal was reported. 

With new experiments and updates underway, the next 
years will become very exciting for axion searches. Both 
helioscope and haloscope experiments are about to probe 
significant parts of complementary parameter regions of 
realistic axion models. For example, an axion signal at 
ADMX would support the hypothesis of i7dm provided by 
axions from the misalignment mechanism. In contrast, an 
axion observation at CAST would point to axion hot / warm 
dark matter which can only provide a minor fraction of 
i7dm so that the dominant contribution can still be pro- 
vided, e.g., by the LSP if SUSY is reahzed in nature. 
In fact, if the axino — the fermionic superpartner of the 
axion — is the LSP, collider signatures predicted for the 
axino LSP [190] could become an additional hint towards 
the existence of the axion and the solution of the strong 
CP problem proposed by Peccei and Quinn. 
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3 Neutralino Dark Matter 

In this section we consider SUSY scenarios in which the 
Hghtest neutrahno Xi is the LSP. This hypothetical par- 
ticle is probably the most studied and most popular SUSY 
dark matter candidate and a concrete example for a weakly 
interacting massive particle (WIMP). Before discussing Xi 
dark matter scenarios, let us review some generic proper- 
ties of SUSY extensions of the Standard Model. For de- 
tails, we refer to dedicated reviews on SUSY [30 l l3T 1[32ll33t 

mm- 

Extending the Standard Model with SUSY, there is a 
superpartner of each Standard Model particle and an ex- 
tended Higgs sector with a least two Higgs doublets. The 
couplings of these superpartners arise by supersymmetriz- 
ing the Standard Model couplings and are thus fixed by 
symmetry. This allows for model independent SUSY pre- 
dictions that can be tested in collider experiments [191|, 
I192iil93 . The masses of the Standard Model superpart- 
ners are governed by the Higgs-higgsino mass parameter 
II and by soft SUSY breaking parameters which depend 
on the SUSY breaking mechanism and thereby on physics 
at high energy scales such as the one of grand unification 
Mqut — 2x 10^^ GeV. The experimental determination of 
the SUSY mass spectrum and of the Higgs masses at col- 
liders can thus provide insights into high scale physics and 
into the SUSY breaking mechanism [194..195,..196.,197J . 

Assuming that SUSY is realized not only as global 
but as a local symmetry [3D], the gravitino G appears 
as the spin-3/2 superpartner of the graviton in addition 
to the Standard Model superpartners. The gravitino is 
the gauge field associated with local SUSY transforma- 
tions and a singlet with respect to the gauge groups of the 
Standard Model. Its interactions — given by the supergrav- 
ity Lagrangian |198[[30] — are suppressed by the (reduced) 
Planck scale [2] 

18, 



Mp 2.4 X 10'" GeV 



Once SUSY is broken, the extremely weak gravitino inter- 
actions are enhanced through the super-Higgs mechanism, 
in particular, at energy/mass scales that are large with 
respect to the gravitino mass m^. Nevertheless, the grav- 
itino can be classified as an extremely weakly interacting 
particle (EWIP). Since the gravitino G is the gauge field of 
local SUSY, its mass is governed by the scale of SUSY 
breaking and can range from the eV scale to scales beyond 
the TeV region [3n[#lll99ll2()()ll2()lll2()2ll2()3ll2()4j . For ex- 
ample, in gauge-mediated SUSY breaking schemes [199[ 
I200p201] . the mass of the gravitino is typically less than 
1 GeV, while in gravity-mediated schemes [5T1I55] it is ex- 
pected to be in the GeV to TeV range. In fact, SUSY 
scenarios in which the gravitino is the stable LSP are well- 
motivated and and will be discussed in the next section. 
In this section, is assumed to be above the neutralino 

mass m^o. This implies an unstable G which can be as- 
sociated with an additional source of Xi dark matter (cf. 
Sect. I3.1| ) and with restrictive BBN constraints on the re- 
heating temperature Tr, (cf- Sect. 13. 2| ). 



The lightest neutralino Xi appears in the minimal su- 
persymmetric Standard Model (MSSM) as the lightest 
mass eigenstate among the four neutralinos being mix- 
tures of the bino B, the wino W, and the neutral higgsi- 
nos -ff" and H^. Accordingly, Xi is a spin 1/2 fermion with 
weak interactions only. Its mass m^o depends on the gaug- 
ino mass parameters Mi and M2 , on the ratio of the two 
MSSM Higgs doublet vacuum expectation values tan/3, 
and on the Higgs-higgsino mass parameter fi. Expecting 
m^o 0(100 GeV), X? is classified as a WIMP. 

Motivated by theories of grand unification and super- 
gravity |205| and by experimental constraints on flavor 
mixing and CP violation [2], one often assumes universal 
soft SUSY breaking parameters at the scale of grand uni- 
fication Mqut; cf. [33ir34.35.2O6] and references therein. 
For example, in the framework of the constrained MSSM 
(CMSSM), the gaugino masses, the scalar masses, and the 
trilinear scalar interactions are assumed to take on the 
respective universal values mi/2, mo, and ^0 at Mqut- 
Specifying mi/2, mo, Ao, tan/?, and the sign of /j,, the 
low-energy mass spectrum is given by the renormalization 
group running from Mqut downwards. 

For ^0 — 0, for example, the lightest Standard Model 
superpartner — or lightest ordinary superpartner (LOSP) — 
is either the hghtest neutrahno Xi or the fighter stau ri, 
whose mass is denoted by mr^ ■ If the LSP is assumed to 
be the LOSP, the parameter region in which m^^ < m^o is 
usually not considered because of severe upper limits on 
the abundance of stable charged particles [2]. However, 
in gravitino/axino LSP scenarios, in which the LOSP is 
the next-to-lightest supersymmetric particle (NLSP), the 
Ti LOSP case is viable and particularly promising for col- 
lider phenomenology as will be discussed in Sects. H] and O 
In Fig. d] (from ^U7\ ) the dotted (blue) fines show 
contours of ttilosp in the (mi/2,TOo) plane for Ao = 0, 
/i > 0, tan/3 = 10. Above (below) the dashed fine, m^o < 
(26) {m^-i < m^o). The medium gray and the light gray 

regions at small mi ^2 are excluded respectively by the 
mass bounds m~± > 94 GeV and toh > 114.4 GeV from 



chargino and Higgs searches at LEP [2]. It can be seen 
that TOjjro = 0(100 GeV) appears naturally within the 

cmssm' 

Before proceeding, let us comment on other poten- 
tial LOSP/LSP/NLSP candidates. For Ao ^ and in 
less constrained frameworks such as models with non- 
universal Higgs masses (NUHM), there are parameter re- 
gions in which the LOSP is the lighter stop h [2TT 1 I2T2 1 
[2T3] or the lightest sneutrino [2T4 l l215 p 216 p 217 p 218 1 
1219] . In fact, since the lightest sneutrino i^i is electri- 



cally neutral and color neutral, the Vi LSP looks at first 
sight like another promising WIMP dark matter candi- 
date within the MSSM. It turns out however that its cou- 
plings (and in particular the one to the Z-boson) are "too 
strong." From the invisible Z-boson width, we know that 
the sneutrino must have a mass mj;^ > Afz/2, where the 
its relic density is typically well below i7dm and/or its in- 
teractions with nuclei are such that it should have already 
been observed in direct dark matter searches (assuming 
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tan/3 = 10, ylo = 0, /i > 




500 1000 
mi/2 [GeV] 

Fig. 4. Contours of mLOSP (dotted blue lines) and Ilosp 
(solid black lines) in the (^1/2,^10) plane for Ao = 0, /j, > 0, 
tan (3 = 10. Above (below) the dashed line, m-o < 
{m:f-^ < m~o). The medium gray and the light gray regions 

show the LEP bounds m-± > 94 GeV and mn > 114.4 GeV, 

^1 

respectively [2]. The contours are obtained with the spectrum 
generator SuSpect 2.34 [208) using rrit = 172.5 GeV and 
TOb(mb)^^ = 4.23 GeV, and with micrOMEGAs 1.3.7 
\2W\ . From [207] . 



a standard dark matter halo profile); cf. Figs. 1 and 2 in 
Ref. [220]. Thus, the MSSM sneutrino LSP is not consid- 
ered a viable dark matter candidate [221j but it may well 
be the NLSP, e.g., in a gravitino/axino LSP scenario [2221 
[^mriEMlETHim^ Fl The lighter stop h is not viable 
as a stable LSP due to severe constraints on exotic sta- 
ble colored particles but another NLSP candidateP^ 
While the NLSP governs cosmological constraints and ex- 
perimental prospects in gravitino/axino dark matter sce- 
narios in a crucial way (cf. Sects. H] and [5]), it can also 
be important in a neutralino dark matter scenario (e.g., 
through coannihilation processes) as will become clear be- 
low. 



3.1 Primordial Origin 

The Xi's were in thermal equilibrium for primordial tem- 
peratures of T > Tf ~ TO^o/20. At Tf, the annihilation 



Variants of the MSSM sneutrinos are actively pursued as 
dark matter candidates; see [220[ and references therein. With 
an admixture of less strongly interacting "right-handed" sneu- 
trinos, sneutrino dark matter interactions can become compat- 
ible with fij;-^ ~ ^2dm and with constraints from direct searches. 

Note that the mass of a long-lived ti NLSP has to re- 
spect the collider bound, m^^ > 250 GeV, inferred from SUSY 
searches at the Fermilab Tevatron 12131. 



mo [GeV] 
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Fig. 5. Regions (black) with 0.087 < O^^'^h'^ < 0.138 in 

the (mi/2,wio) plane for Aq — 0, ^ > 0, tan/9 — 10, and mt — 
172.7 GeV. In the dark gray triangular region, m~o > irir^ . The 
light gray region at small mi/2 is excluded by the requirement 
of correct electroweak symmetry breaking or by sparticle search 
limits [228j , the two medium shaded (light pink) bands by the 
LEP bound mn > 114 GeV, and the small light shaded (green) 
spot by the 6 — > 57 constraint: 2.65 < BR(fe S7)/10^'' < 
4.45. The dark shaded (red) band is compatible with having a 
Standard-Model-like Higgs boson near 115 GeV. From [228[ . 



rate of the (by then) non-relativistic xT^ becomes smaller 
than the Hubble rate so that they decouple from the ther- 
mal plasma. Thus, for T < T{, their yield y^ro e 



is given by Fiff 



dec 



n^o / s 

A.1 ' 

where n^~if^ is the (equi- 
librium) number density of x?'s. Depending on details 
of the Xi decoupling, Y~o'^ is very sensitive to the mass 

spectrum and the couplings of the superparticles. Indeed, 
convenient computer programs such as DarkSUSY |225j or 
micrOMEGAs |2091I2101I2251[^ are available which allow 
for a numerical calculation of LOSP decoupling and of 
the resulting thermal relic abundance in a given SUSY 
model. 

The l^Logp contours shown by the solid black lines in 
Fig. H] illustrate that the Xi LSP yield can easily vary by 
more than an order of magnitude. Because of this sensi- 
tivity, the associated thermal relic density 



(therm 1 2 



m^oY^fs{To)hyp, 



(27) 



agrees with only in narrow regions in the parameter 

space; pc/[s(To)h^] = 3.6 x lO^^ GeV [Ij. This can be seen 
in Fig. [5] (from f228]) where the black strips indicate the 
region with 0.087 < l2*JS°™ft2 < q ^sg. 

Remarkably, it is exactly the small width of the re- 
gions with J7io''''™ = fidm, which could help us to identify 
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Xi dark matter. If sparticles exist and if they are produced 
at colliders, the data analysis will aim at determining the 
SUSY model realized in nature 196,197J. For the recon- 

thorm possible 

assuming a standard thermal history of the Universe. Be- 
cause of the sensitivity of !7io'^'''° with respect to the SUSY 

model, an agreement of the obtained fi^g 



structed model, a precise calculation of i7|p' 



therm 



with i7dm will 

then be a strong hint for the Xi LSP providing J7dm and 
for a standard thermal history up to the Xi"decoupling 
temperature Tf. Since Xi's decouple already as a non-re- 
lativistic species, it is also guaranteed that they are suffi- 
ciently cold to allow for cosmic structure formation. 

In fact, f2^°™ exceeds f^dm in most of the parame- 

ter space. Thus, regions with Sl^^""^ ~ fidm are some- 

what special and associated with particularly efhcient neu- 
tralino annihilation in the early Universe. Depending on 
the origin for the efficient annihilation, these regions can 
be classified as follows; cf. [MllHSl l^l^l^ and refer- 
ences therein: 

— Bulk region: This region is associated with light slep- 
tons li, mj^ < 200 GeV, so that neutralinos can an- 
nihilate efficiently via slepton exchange into a lepton 
pair: XiXi ■ This region is often in tension with 
the LEP Higgs bound. For example, in the CMSSM 
scenarios considered in Figs. [4] and [5l the bulk region 
appears around (mo,mi/2) — (60 GeV, 200 GeV). 

— Focus point region/hyperbolic branch: This region is 
associated with a Xi with a significant higgsino admix- 
ture so that 

Xlxi ^ W+W-, Z°Z° become efficient. 
(In fact, for a purely bino-like neutralino, Xi = 
annihilation into these final states cannot occur.) In 
Fig.O the black region with /2io'"^™ — ^dm at large ttiq 
is associated with these neutralino annihilation chan- 
nels being efficient. 

— Coannihilation region: In this region, the NLSP has 
a mass very close to the one of the neutralino LSP. 
Thereby, the number density of the NLSP during the 
freeze out of the Xi LSP is sizable and Xi~NLSP coan- 
nihilation processes can enhance the efficiency of neu- 
tralino annihilation. In Fig. [5l the black region with 
l^thcrm ^ -^^ above the dark gray stau LOSP 

Xi 

region is associated with efficient Xi-'?L-coannihilation 
processes. Moreover, in regions in which m^_^ is close to 
m^o, also Y?~ti coannihilation can allow for i7lo°''™ ~ 

^dm ■ 

— Higgs funnel: In this region, 2 rn^o is close to the mass 
of the CP odd Higgs boson to^o 2m^a, so that 
neutralino annihilation proceeds very efficiently via the 
A° resonance. For large tan/3, the following annihila- 
tion channel becomes particularly efficient due to an 
tan /3-enhanced coupling of the b quark: XiXi ~^ 
A° bb. 



In addition to /?io 

Xl 



:)thcrm fj-Qjjj thermal freeze out, Xi's can 

also be produced non-thermally in late decays of graviti- 
nos. Because of their extremely weak interactions, unsta- 



Table 1. Assignments of the index i, the gauge coupling gi, 
and the gaugino mass parameter Mi, to the gauge groups 
U(1)y, SU(2)l, and SU(3)c, and the constants ki, yi, and uji. 



gauge group 


i gi 


M, 


ki 


(y./io-^^) 




U(1)y 


1 9' 


Ml 


1.266 


0.653 


0.018 


SU(2)l 


2 9 


Ma 


1.312 


1.604 


0.044 


SU(3)c 


3 gs 


Ms 


1.271 


4.276 


0.117 



ble gravitinos with m^o < mg < 5 TeV have long life- 
times, Tg > 100 s (cf. Fig. 1 of Ref. [63J), and decay typi- 
cally during or after BBN into the LSP and into Standard 
Model particles. While the decay into the LSP can pro- 
ceed either directly or via a cascade, each gravitino decays 
into one LSP. Thus, the resulting non-thermally produced 
(NTP) neutralino density is given by 



NTP J 2 



h m^o Yg s(To) h / pc , 



(28) 



where Yg = ng/s denotes the gravitino yield prior to de- 
cay. Although gravitinos with mg > m^o = 0(100 GeV) 
are extremely weakly interacting and not in thermal equi- 
librium with the primordial plasma, they can be produced 
efhciently in thermal scattering of particles in the hot 
plasma|i3 Derived in a gauge-invariant treatment, the re- 
sulting thermally produced (TP) gravitino yield at a tem- 
perature Tiow < Tr reads p3nr232, 207, 235] 



(-^low) 



i=l 



X In 



.9i(7r 



3m| J 
1010 GeV 



, (29) 



with Ui, the gauge couplings gi, the gaugino mass param- 
eters Mi, and ki as given in Table [TJ Here Mi and gi 
are understood to be evaluated at the reheating temper- 
ature after inflation Tr 232 Using = ^^^(Tlow) 
as given by ((29|) . one finds that H-i^^ is sensitive to the 

Xl 

gravitino mass mg and to the reheating temperature Tr. 
Thus, the dark matter constraint n^Jl""'' + H^^^ < Hdm 

' Xl Xl — "™ 

does imply an upper limit on Tr. This limit is particularly 
restrictive in scenarios with i7lo°''™ ~ J7dm [235 and/or 

TO- <C M^(Tr). Since non-thermally produced x^'s can be 
hot/warm dark matter, additional constraints from LSS 
and potential solutions to small scale structure problems 



can occur, in particular, for f2~o 



In this review I do not discuss gravitino production from 
inflaton decays which can be substantial depending on the in- 
flation model; see, e.g., [229ii230j . 

Note that the field-theoretical methods applied in the 
derivation of ^ [23T][232][207l[233] require weak couplings 
gi <^ 1 and thus T ^ 10® GeV. For an alternative approach, 
see [2341. 



NTP 
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Let us comment at this point on the Tr definition as 
discussed in Ref. [20 7j . The analytic expression (f29|) is de- 
rived by assuming a radiation-dominated epoch with an 
initial temperature of Tr |23ip232j . In a numerical treat- 
ment, the epoch in which the coherent oscillations of the 
inflaton field dominate the energy density of the Universe 
can also be taken into account, where one usually defines 
Tr in terms of the decay width of the inflaton field [5T1 
12071 . In fact, the numerical result for ^^^(riow) agrees 
with the analytic expression (f29|) for |207) 



90 



1/4 



1.8 



(30) 



which satisfies ~ 1.8ifrad(7R) with the Hubble pa- 
rameter i/rad(r) = v/g,(T)7rV90rVMp and an effec- 
tive number of relativistic degrees of freedom of 5, (Tr ) = 
228.75. Thus, ^ provides the Tr definition to which 

the yield ([29|l applies. For an alternative definition T^^' = 
[90/ (5* (4*') )] 1/4 yj^MFTJ given by = ei?rad(T|l), 
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Fig. 6. Upper limits on the reheating temperature Tr af- 
ter inflation as a function of the gravitino mass mg in the 



the associated numerically obtained yi^(riow) is repro- CMSSM with xi dark matter for 



m 



duced by the analytical expression 

Tr with ./UUT^^ 



after substituting 



1/2 



300 GeV, mo 



[207], 



Figure [6] (from [63j) shows the Tr limit imposed by 
NTP/i2 < 0.118 (dotted line labeled as "i7LSp") asafunc- 



xi 

tion of for the CMSSM scenario with mi/2 = 300 GeV, 
mo = 141 GeV, Aq = 0, and tan/3 = 30. Here the Tr 
limits apply to the Tr definition associated with ^ = 3, 
i.e., T^ ~ 3iJrad(2R,). In the shaded (light-orange) re- 
gion, TOg < m^a — 117 GeV. The Tr limit becomes 
more restrictive for mg — *■ m^o due to the mg-dependent 
goldstino component in (|29p . For large mg, this spin- 
1/2 component becomes negligible and the Tr limit is 
governed by the m^-independent spin-3/2 contribution 
in Note that the shown limit is conservative since 

the thermal relic density associated with this SUSY sce- 
nario, n^^°™h^ ~ 0.111, is not taken into account. 



3.2 Cosmological Constraints 



141 GeV, j4o = 0, tan /3 = 30 and for the Tr definition as- 
sociated with ~ 3-?/rad(rR). The associated Xi LSP mass 
is m-o = 117 GeV — as indicated by the shaded (light-orange) 
region in which mg < m~o — and the thermal relic density 



thcrm^2 



0.111. Above the dotted line labeled as "i^LSp" 



Late decaying gravitinos are not only associated with the 
contribution to f^^oh^ but also with the injection of 
energetic Standard Model particles. Because of the long 
G lifetime of > 1 s for mg < 20 TeV [63], those decay 
products are emitted during/after BBN and can thus af- 
fect the abundances of the primordial light elements [52) 
[MUSOMIIIMSIISIISS]- In fact, this is a concrete (and prob- 
ably the most prominent) example for the non-thermal 
BBN-affecting processes mentioned in the Introduction. 

The dominant mechanism affecting BBN depends on 
Tq, or, more generally, on the time t at which the electro- 
magnetic or hadronic energy is injected into the Universe. 
For Is < t < 100 s, energetic hadrons are stopped effi- 
ciently through electromagnetic interactions so that the 
direct destruction of light elements is subdominant. The 
presence of additional slow hadrons still can change the 



the Xi density from decays of thermally produced gravitinos 

exceeds H^S^^h^ = 0.118. The lines labeled as D, ^He, "He, ^Li, 
^1 

and '^Li are upper limits on Tr (95% CL) inferred from obser- 
vationally inferred primordial abundances of the corresponding 
light elements. Above those limits, BBN is reprocessed in an 
intolerable way by the Standard Model particles emitted in 
late decays of thermally produced gravitinos. From [63] . 

ratio of protons to neutrons through interconversion pro- 
cesses and thus affect the abundance of the light elements. 
For 100 s < t < 10''s, energetic hadrons and, in particu- 
lar, neutrons cannot be slowed down significantly. Accord- 
ingly, they can reprocess efficiently the produced light el- 
ements through hadrodissociation processes. The effect of 
electromagnetic energy release is negligible for t < 10^ s 
as the interaction with the background particles thermal- 
izes quickly any high-energy photons or leptons emitted in 
the gravitino decay. Towards later times, electromagnetic 
energy release becomes important. For 10^ s < i < 10^^ s, 
the reprocessing of light elements through energetic elec- 
tromagnetic showers, i.e., photodissociation, can become 
more significant than hadrodissociation. (For more details, 
see, e.g., [541161] and references therein.) 

Including these mechanisms in calculations of BBN, 
observationally inferred abundances of primordial D, '^He, 
^He, ^Li, and ''Li have been used to provide limits on 
quantities such as [SWCTIIS^ 

^cin,had = ^cm,had 7 (^1) 

where ecm.had is the (average) electromagnetic/hadronic 
energy emitted in a single G decay. While the yield prior 
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to decay is given for thermally produced gravitinos 
by (I29p . eem.had depends strongly on the sparticle spec- 
trum. Once eom.had is calculated for a given SUSY model, 
^em,had limits Can basically be translated into m ^-de pen- 
dent upper limits on the reheating temperature |2351IS5] . 

For the exemplary CMSSM point TOi/2 = 300 GeV, 
mo — 141 GeV, Aq — 0, and tan/3 = 30, the mg-dependent 
BBN constraints on (95% CL) are shown in Fig. [6] 
(from [B3J) for the Tr definition given by ~ 3i/i.ad(7R). 
The curves are inferred from observationally inferred abun- 
dances of D, '^He, ''He, ^Li, and ''Li (as labeled). The 
^He limit governs the BBN constraints for > 7 TeV 
{tq < 10^ s) where proton-neutron interconversion can af- 
fect BBN. The D ('^He) limit is the most restrictive one in 
the region in which the constraints from hadrodissociation 
(photodissociation) are most relevant, 0.3 TeV < < 
7 TeV (mg < 0.3 TeV) or 10^ s < rg < 10^ s (rg > 10^ s). 
For > 40 TeV, the BBN bounds disappear since the 

G's decay well before the onset of BBN, i.e., <C 1 s. 

The range of allowed values of the reheating temper- 
ature is crucial for our understanding of inflation and 
for the viability of potential explanations of the matter- 
antimatter asymmetry in our Universe. For example, ther- 
mal leptogenesis with hierarchical heavy right-handed Ma- 
jorana neutrinos — which provides an attractive explana- 
tion of this asymmetry — requires very high reheating tem- 
peratures of Tr > 10^ GeV [238[|239y240y24iy242j and 
thus > 7 TeV for the SUSY model considered in 
Fig. [51 For smaller tUq, the Tr limit can be as restrictive 
as Tr < 10^ GeV which is known as the gravitino problem. 
Note that flavor effects [243l[244l[24Tl[2i2] do not change 
the lower bound Tr > 10^ GeV required by successful 
thermal leptogenesis with hierarchical right-handed neu- 
trinos [24ip242j . However, there are special settings [2451 
I246P247] that allow for a CP asymmetry above the Davidson- 
Ibarra bound |239| and for a relaxed lower Tr bound. 
Moreover, for (nearly) mass-degenerate heavy right-handed 
Major ana neutrinos, resonant leptogenesis can explain the 
matter-antimatter asymmetry at smaller values of Tr [2481 
I249[l250|l251j . Another example for a framework in which 
the limit Tr > 10^ GeV is relaxed is non-thermal leptoge- 
nesis; see, e.g., [252] and references therein. 

Before proceeding one should stress that the Tr lim- 
its inferred from f2~^^ < fidm and from BBN rely cru- 

daily on assumptions on the cosmological history and the 
evolution of physical parameters. For example, for a non- 
standard thermal history with late-time entropy produc- 
tion, the thermally produced gravitino yield can be diluted 
^ Yg/S by a factor 6 > 1 [203 so that T^^'^ ^ ,5 T^^"". 
Moreover, T^'^^ can be relaxed if, e.g., the strong cou- 
pling gs levels off in a non-standard way at high tempera- 
tures |253j . This emphasizes that the Tr limits discussed 
above rely on the assumptions of a standard cosmological 
history and of gauge couplings that behave at high tem- 
peratures as described by the renormalization group equa- 
tion in the MSSM. While tests of these assumptions seem 
inaccessible to terrestrial accelerator experiments, the fu- 
turistic space-based gravitational-wave detectors BBO or 



DECIGO [87] — mentioned already in the Introduction — 
could allow for tests of the thermal history after inflation 
and could even probe Tr 88,89J in a complementary way. 

3.3 Experimental Searches and Prospects 

For experimental tests of the Xi dark matter hypothe- 
sis, three complementary techniques exist: indirect, direct, 
and collider searches. While there is an enormous activ- 
ity in each of those fields, I will summarize only the main 
ideas. For more detailed discussions, see |254[|255ll2561l257l 
258J and references therein. 

Let us first turn to indirect searches. Since dark mat- 
ter clumps, one expects regions with an increased x? den- 
sity such as galaxy halos, the center of galaxies, and the 
center of stars. While Xi pa-h annihilation well after x^ 
decoupling is basically negligible for calculations of i^^o, 
it should occur at a significant rate in these regions. The 
resulting Standard Model particles should then lead to en- 
ergetic cosmic rays and thereby to an excess of photons, 
neutrinos, positrons, and antiprotons over backgrounds 
expected from standard cosmic ray models without dark 
matter annihilation. In fact, for example, data from the 
Energetic Gamma Ray Experiment Telescope (EGRET) 
has already been interpreted as evidence for Xi annihi- 
lation [259] within SUSY models that will be testable 
in direct and collider searches. For a discussion of these 
and other potential hints, see [257,255,258 and references 
therein. 

In direct searches, one looks for signals of x^'s — or 
more generally WIMPs — passing through Earth that scat- 
ter elastically off nuclei. Being located in environments 
deep underground that are well shielded against unwanted 
background, an enormous sensitivity has been reached by 
a number of experiments [260l[26Tll262l[263l[264l[265l[266l 
1267[l268j . Since no unambiguous signal of a Xi~iiucleus 
scattering event has been observed so far, m^o-dependent 

upper limits on the respective Xi cross section are ob- 
tained. 

Figure [7] shows current limits on the spin- independent 
neutralino-nucleon cross section provided by CDMS [2661 
(red), XENONIO [265] (dark blue), CRESST II [267,268 
(light blue), ZEPLIN II ;264] (dark green), WARP [ 263] 
(medium green), and EDELWEISS I [IHI] (hght green) 
(from bottom to top at a WIMP mass of 80 GeV, as la- 
beled). The light (yellow) patch indicates the region in 
which the DAMA experiment is reporting the observa- 
tion of a signal with the expected annual modulation [2691 
1270] . As can be seen, this DAMA signal region is in a re- 
gion in which null events were observed by ZEPLIN II, 
CRESST II, CDMS, and XENONIO. The DAMA signals 
might thus be interpreted as signals of a "non-standard" 
dark matter candidate (see e.g. [271] ). where cosmological 
and astrophysical constraints can allow for crucial viabil- 
ity tests (see e.g. ^7^); see also |?751[?7il[?75] and refer- 
ences therein. The various patches, crosses, and contours 
within the light gray region labeled as "supersymmetric 
models" indicate the (favored) parameter space of SUSY 
models with a x? LSP [?mi277y^78y^79y^80] . As can be 
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Fig. 7. Limits on the spin-independent neutralino-nucleon 
cross section from direct searches as a function of the neu- 
tralino mass (or WIMP mass). The shown hmits are pro- 
vided by CDMS [266 (red), XENON 10 [263 (dark blue), 
CRESST II [257112681 ' (light blue), ZEPLIN II [2M] (dark 
green), WARP [263] (medium green), and EDELWEISS I [26T] 
(light green) (from bottom to top at a WIMP mass of 80 GeV, 
as labeled). The DAMA signal region [25M270] is indicated 
by the light-shaded patch (yellow, as labeled) and the (fa- 
vored) parameter space of SUSY models with a Xi LSP [7751 
[2771I278II279II280] by the shadings/contours/crosses at WIMP 
masses m^o > 40 GeV (gray/purple/red, as labeled). Stan- 
dard galactic halo parameter are assumed, i.e., a local halo 
density of dark matter of po — 0.3 GeV/cm'^ and a char- 
acteristic halo velocity of vq = 220 — 240 km/s. The fig- 
ure was build by using the dark matter plotter available at 
[http : //dmtools .berkeley . edu/ 1 imitplot s / ^ and maintained 
by R. Gaitskell and J. Filippini. 



seen, the current best limits given by the CDMS [262il266j 
and the XENONIO [265j experiments disfavor already a 
sizable part of the SUSY parameter space; see, for ex- 
ample, [2551l256|l206j and references therein. These limits, 
however, depend on the assumed Xi flux at the detector lo- 
cation. Standard galactic halo parameter are assumed, i.e., 
a local halo density of dark matter of po = 0.3 GeV/cm'^ 
and a characteristic halo velocity of vq — 220 — 240 km/s. 
Indeed, those assumptions are subject to significant uncer- 
tainties due to possible inhomogeneities in the dark mat- 
ter distribution in galaxies. Such inhomogeneities should 
manifest themselves also in indirect searches which can 
help to reduce those uncertainties. Once Xi events are ob- 
served in direct searches, one can succeed in reconstruct- 
ing the Xi velocity distribution [281j . By analyzing the 
recoil spectra, m^o can even be estimated in a way that 
is independent of the dark matter density on Earth [282 . 

In most searches for SUSY at colliders, it is assumed 
that R-parity is conserved. Accordingly, one expects that 



superpartners are produced in pairs before decaying via 
cascades into the LSP and energetic fermions. As a weakly- 
interacting particle, every Xi LSP produced will escape 
the detector without leaving a track. Thus, the existence 
of SUSY and the Xi LSP has to be inferred from studies 
of missing transverse energy E^^^^ and of energetic jets 
and leptons emitted along the cascades. Along these lines, 
ongoing investigations are pursued based on data from pp 
collisions with a center-of-mass energy of = 2 TeV 
at the Fermilab Tevatron Collider. While lower limits on 
the masses of squarks and gluinos have been extracted, 
no evidence for SUSY or the Xi LSP has been reported 
so far |2831l284j . With the first pp collisions with y/s = 
14 TeV at the CERN Large Hadron Collider (LHC) ex- 
pected in the year 2009, there are high hopes that the 
new energy range will allow for a copious production of 
superpartners. Here large E^^^^ will be the key quantity 
for early SUSY searches 



[285[I286| . Despite an enormous 
potential for mass and spin measurements of SUSY parti- 
cles at the LHC [287) . additional precision studies at the 
planned International Linear Collider (ILC) [288p289j ap- 
pear to be crucial for the identification of the x? LSP [2541 
[2901. 



4 Gravitino Dark Matter 

The gravitino G has already been introduced in the previ- 
ous section and its appearance is an unavoidable implica- 
tion of SUSY theories including gravity. In this section we 
consider the possibility of the gravitino LSP which is well 
motivated, for example, in gauge-mediated and gravity- 
mediated SUSY breaking schemes [199,l200l[20T]. Indeed, 
without consensus on the SUSY breaking mechanism and 
the SUSY breaking scale, one may well consider the grav- 
itino mass mg as a free parameter to be constrained by 
cosmological considerations (cf. Sects. [4Tl and [42)1 and by 
collider experiments (cf. Sect. I4.3|) F^ 

Being a singlet with respect to the gauge groups of the 
Standard Model, the gravitino LSP is a promising dark 
matter candidate that can be classified as an EWIP as 
mentioned above. In fact, it must not be massive since 
even a light gravitino (e.g., = 1 keV) can evade its pro- 
duction at colliders because of its tiny interaction strength. 

Let us recall important characteristics of gravitino in- 
teractions: (i) Gravitino interactions are suppressed by in- 
verse powers of Mp. Indeed, for example, gravitino-gau- 
gino-gauge boson couplings are described by dimension 
five operators and an energy scale appears in the numer- 
ator of the respective vertex. Gravitino interactions can 
thereby be enhanced at very high energies, (ii) Through 
the super-Higgs mechanism, the interactions of the spin- 
1/2 components of the gravitino (i.e., the interactions of 
the goldstino components) are enhanced at energy/mass 
scales that are large with respect to the gravitino mass 



In scenarios with R-parity violation, the gravitino LSP can 
decay into Standard Model particles. Then, one may be able 
to infer its mass from the decay spectra possibly observed 
in indirect dark matter searches 12911. 
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mg, i.e., a light gravitino interacts more strongly than a 
heavy gravitino 13 

Considering the case of the G LSP, in which the LOSP 
is the unstable NLSP that decays eventually into the G 
LSP, both cases m^o < and m^^ < m^o (cf. Fig. [4]) 
are viable as already mentioned. In less constrained frame- 
works such as models with non-universal Higgs masses 
(NUHM), also other LOSP/NLSP candidates are still pos- 
sible such as the lighter stop h ^HWIH^ or the light- 
est sneutrino ?i p23[|217il224il218il219j . However, the fi 
NLSP case is probably the most promising one from the 
phenomenological point view and is discussed more exten- 
sively than the other NLSP cases in this review. 



4.1 Primordial Origin 

The potential primordial origin of gravitino dark mat- 
ter depends on the mass that governs its interaction 
strength, on the SUSY model, on the cosmological his- 
tory, on the inflation model, and on the reheating tem- 
perature after inflation. The gravitino LSP can be a ther- 
mal relic j292| or be produced in thermal scattering of 
particles in the primordial plasma [293 1I294U231U232U2071 
[234, . Additional more model dependent gravitino sources 
are NLSP decays [295,296,297,223] and decays of scalar 
fields such as the inflaton ^229,i230]. The latter produc- 
tion mechanism is not discussed in this review but can be 
substantial depending on the inflation model. 

Light gravitinos can have sufficiently strong interac- 
tions for being in thermal equilibrium with the primor- 
dial plasma. For example, for < 2 keV, the gravitino 
decoupling temperature is Tf < 1 TeV; cf. [34]. Before 
gravitinos decouple as a relativistic species at Tf, i.e., for 
Tr > T > T{, the spin-1/2 components of the gravitino 
were in thermal equilibrium. Thus, their "hot" thermal 
relic density is the one of a spin-1/2 Majorana fermion 



therm L 2 



0.115 



100 

g*s{T{) 



- 100 eV 



(32) 



and the right amount of dark matter, if^^^^ ~ ^2dm, is 
provided for ~ 100 eV and g*s(Tf) ~ 100. However, 
the present root mean square velocity of such gravitinos 



/ rms,0^ n -7-7 ^"""^ I 



thcrmj,2\ evy/^ 



exceeds significantly the constraints from observations and 
simulations of cosmic structures listed, e.g., in Table 1 of 



Ref. [298]. In fact, for l?^^"" 



i7dm, these constraints im- 



ply niQ > 500 eV and g*siTf) > 500 which is well above 
the 228.75 degrees of freedom of the MSSM. Indeed, for 



Note that this interaction-strength dependence on the mass 
is different than in the axion case, i.e., a light axion is less 
strongly interacting than a heavy axion; cf. Sect. (2] 

A long-lived stop ti NLSP is not feasible in the 
CMSSM [2T21I2T3]. 



9*s{Ti) = 228.75, this scenario is excluded by the dark 
matter constraint < f2dm once a standard cosmologi- 
cal history is assumed. With a non-standard thermal his- 
tory, light gravitinos can still be viable thermal relics if 
their abundance is diluted by entropy production, which 
can result, for example, from decays of messenger fields 



in gauge-mediated SUSY breaking scenarios |299P300P30H 
1^)2.303.304] . 

Gravitinos with mg > 0.1 GeV (1 GeV) have a high 
decoupling temperature of Tf > 10" GeV (lO^^ GeV) [207] 
because of their extremely weak interactions. Thus, those 
gravitinos have never been in thermal equilibrium with 
the primordial plasma even for a reheating temperature 
as high as Tr ~ 10^" GeV. At high temperatures, how- 
ever, those gravitinos can be produced efficiently in ther- 
mal scattering of particles in the primordial plasma, as al- 
ready discussed for the case of an unstable G in Sect. 13.11 
In the case of the stable G LSP, the thermally produced 
gravitino yield ^^^(To) = Y^^iTio^) given in ((29|) leads 
to the following thermally produced (TP) gravitino den- 
sity [23ni232ll207j 



nfh^ = r|P(To) s(To) h^/p. 



3m% 

Cr 



In 



( \ 

VlOO GeV/ 



100 GeV/ V IQiOGeV 



Tr 



(34) 



with uji as given in Table [TJ For the case of universal gaug- 
ino masses Mi^2,3 — ^■1/2 at Mqut and ^ Mi, i.e., 
(1 + Mfl 2,171%) ~ Mf/3m%, n%^h'^ can be approximated 

Cr Cr Cr — — 

by the convenient expression [74]^ 



0.32 



^ 10 GeV N,/ mi/2 in ^ 

y^^JKTTev) Vl08 GeV J- ^^^^ 



The thermally produced gravitinos do not affect the ther- 
mal evolution of the LOSP (or NLSP) prior to its decay 
which occurs typically after decoupling from the thermal 
plasma. Moreover, since each NLSP decays into one G 
LSP, the NLSP decay leads to a non-thermally produced 
(NTP) gravitino density [295ll296ll306ll223] 



'"g ^NLSP 



s{To)h^/pc 



(33) so that the guaranteed gravitino density is given by 



(36) 



(37) 



Here the one-loop evolution described by the renormaliza- 
tion group equation in the MSSM is used to evaluate gi{T^) at 
a representative scale of Tr = 10* GeV and to express M^Tr) 
at that scale in terms of Mi,2,3(Mgut) = '7ii/2- Note that 
the use of the two-loop evolution will lead to a somewhat 
smaller prefactor. Going in the other direction, i.e., using the 
two- loop evolution to express M3(rR) in terms of the physical 
gluino mass, a prefactor has been found that is about twice as 
large [305] as the one obtained when the one-loop evolution is 
used to express M3(Tr) in terms of M3(l TeV) [23l] . 
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While is sensitive to Mi and Tr for a given mp,, 

i?-"^^ depends on i^NLSP ~ ^losp ^'^'^ thereby on details 
of the SUSY model realized in nature; cf. Sect. 13.11 For 
the case of a charged slepton li — such as the lighter stau 
Ti — being the NLSpIH simple approximations have been 
used such as pMll^l^lgggi lTlirfg] 

Ijf ^ 0.7 X 10-1^ (^) , (38) 

where Yj = nr Is and n? denotes the total U number 
density assuming an equal number density of positively 
and negatively charged Zi's. Note that the yield ([38| is 
in good agreement with the curve in Fig. 1 of Ref. 296j 
that has been derived for the case of a purely 'right- 
handed' Ti ~ Tr NLSP with a mass that is significantly 
below the masses of the lighter selectron and the lighter 
smuon, ^ rrie-^^-p-^, and with a bino-like lightest neu- 
tralino, Xi — that has a mass of mj^ = 1.1 rrirj^. In the 
case of an approximate slepton mass degeneracy, to?j < 
'^ei,/ii ^ l.lm^Tj, the Ti NLSP yield (|38|) can become 
twice as large due to slepton coannihilation processes [2961 
1207] . Approaching the x?""?! coannihilation region, m^o « 
mfj, even larger enhancement factors occur; cf. Fig.[4l On 
the other hand, while existing studies of Yj^'^'^ focus mainly 

on the h Ir case [296 l l222 p 207 | [308 ] . it has recently been 
found that a sizable left-right mixing of the stau NLSP 
can be associated with an increase of its MSSM couplings 
and thus with a significant reduction of l^*"^ [31011311] . 
The study [311] shows also explicitly that stau annihila- 
tion at the resonance of the heavy CP even Higgs " can 
be particularly efficient via tiTi — > bb and thus be 

associated with exceptionally small values of Y-^'^'^. 

Figure E (from [SJ) shows that r~^<='= < 4 x lO^^^ 
occurs in special regions even within the CMSSM. The 
points B and C mark the regions in which Y~°'^ is excep- 
tionally small due to stau annihilation at the reso- 
nance and due to enhanced stau-Higgs couplings leading 
to efficient annihilation into Higgs bosons, respectively. 
Within these regions, i7~"^^ is negligible and otherwise 
restrictive cosmological constraints can be evaded as will 
be explained in more detail in the next section. 

An exceptional reduction of can occur also in a 

11 

non-standard thermal history with late-time entropy pro- 
duction after the decoupling of the li NLSP and before 
BBN [T13. 207.68J or in low Tr scenarios [TT]- Focussing 
on a standard cosmological history, we disregard such pos- 
sibilities and consider in the following mainly the more 
generic Y~^°'^ values described approximately by 

Let us proceed with the discussion of the relic grav- 
itino density. Indeed, scenarios with fi^ — J7dm are found 
for natural mass spectra and for a wide range of m^-T ^^ 
combinations. This is illustrated in Figs.|9l[T0l andllT] 

For yN£sp in the sneutrino and stop NLSP cases, see 
Refs. [222]|223y220ll2T9] and [222]|307II212H308] . respectively. 
Reference [308) covers also the stau NLSP case for ri ~ tr 
and the effects of Sommerfeld enhancement on Vnlsp- 



tan/9 = 55, Ag = 2mo, fi > 




400 1000 2000 4000 



mi/2 [GeV] 

Fig. 8. Contours of Y^^'' (as labeled) in the (mi/2, rno) plane 
for tan/3 = 55, Aq = 2mo, and n > 0. Darker shadings imply 
smaller Y-^'^ values. The dashed lines are contours of m^pj^ — 
100, 300, and 600 GeV (from left to right). The large light- 
shaded region in the lower left corner is excluded by bounds 
from direct Higgs and SUSY searches (or by the appearance of 
a tachyonic spectrum) . In the region to the left of the vertical 
solid and dotted lines, m^o < 114.4 GeV [2] and B{b S7) > 
4.84 X 10~* [312] . respectively. In the white area, m~o < m^^^. 
From [3TT] . 

In G LSP scenarios, upper limits on Tr can be de- 
rived since < l?dm [293 p 296 p 314 p 298 p 207 p 315 l 3T6] . 
These wg-dependent limits are shown in Fig. [9| (from 
[207] ) for the Tr definition dSOj). They can be confronted 
with inflation models and with potential explanations of 
the cosmic matter-antimatter asymmetry in the same way 
as in the case of the Xi LSP (cf. Fig. O in Sect. [3l2 
For a given nj,^ and a h NLSP with jSHj), the bound 
I^NTP ^ fidm — gives upper limits on mg and mj^. 
Those limits are shown by the thin solid lines in Fig. [TOl 
(from [72]). In Fig.[n](from 'Wl]) regions with % e H^'^ 
are shown for Tr = 10^, 10^, and 10^ GeV, where Tr is 
defined as in ([30ll . Here both J7^^ and are taken 
into account for — mo within the CMSSM. 

While thermally produced gravitinos have a negligible 
free-streaming velocity today, gravitinos from NLSP de- 
cays can be warm/hot dark matter. In the ti NLSP case, 
for example, upper limits on the free-streaming velocity 
from simulations and observations of cosmic structures ex- 
clude nif^ < 0.7 TcV for I^^^f^ [gggj. Such scenarios 

(gray band in Fig. [TO]), however, reauirj^ m^^ > 0.7 TeV 



Note that the BBN constraints discussed below can point 
to a Ti NLSP mass of ^ > 2 TeV for J7^'^^~ J?dm as shown 
in Fig. 1101 In the part of the gray band in which mf ^ < 3 TeV 
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Fig. 9. Upper limits on the reheating temperature Tr in the 
G LSP case, which are associated with the Tr definition (|30|) . 
On the upper (lower) gray band, J?g^ € for Mi, 2, 3 = 

mi/2 = 500 GeV (2 TeV) at Mgut- The corresponding limits 
from fl^h^ < 0.126 shown by the dashed and dotted lines are 
obtained respectively with ((M]) for Mi/ 10 = M2/2 = M3 = 



mi/2 at Mgut and with the result of Ref. [231| for M3 = mi/2 
at Mgut- From [207] , 



anyhow and could even resolve the small scale structure 
problems inherent to cold dark matter |318[l319il320j . 



4.2 Cosmological Constraints 

In the G LSP case with conserved R-parity, the NLSP 
can have a long lifetime tnlspH This is illustrated by the 
dotted tnlsp contours in Figs. [10] and [TlJ In particular, 
for a li NLSP, one finds in the limit m; ^ 0, 




10 100 
[GeV] 



1000 



Fig. 10. Cosmological constraints on the masses of the grav- 
itino LSP and a purely 'right-handed' h NLSP with pg)) . 



The gray band indicates fi 



NTP 



Above this band. 



> 0.126. On the thin solid lines labeled with / val- 
ues only / i7dm is provided by ff~^^ . The dotted lines show 
contours of Tf_^. Due to CBBN, the region below the solid and 
the long-dash-dotted (red) lines is disfavored by observation- 
ally inferred abundances of ^Be and ®Li, respectively [79] . The 
effect of electromagnetic and hadronic energy injection on pri- 
mordial D disfavors the regions inside the short-dash-dotted 
(blue) curves and to the right or inside of the short-dashed 
(blue) curves, respectively. Those curves are obtained from the 
severe and conservative upper limits defined in Sect. 4.1 of [298| 
based on results of [591161) . The region below the dashed (green) 
line is disfavored by the effect of electromagnetic energy injec- 
tion on "^He/D FSTI. While the constraints from hadronic en- 
ergy injection are obtained for a purely 'right-handed' li ~ Zr 
NLSP, the ones from electromagnetic energy injection are valid 
for the Ti NLSP case with a visible electromagnetic energy of 
Evis ~ eom = O.SEt released in ri Gt. From [79] . 



r-\ii ^ Gl) = 



487rTO|M|, 




^ , (39) 



which holds not only for a charged slepton NLSP but also 
for the sncutrino NLSP vi. Expressions for the lifetimes of 
the Xi NLSP and the ti NLSP are given, e.g., in Sect. IIC 
of Ref. [223] and in Sect. 2.2 of Ref. [212], respectively. 

If the NLSP decays into the G LSP occur during or af- 
ter BBN, the Standard Model particles emitted in addition 
to the gravitino can affect the abundances of the primor- 
dial light elements as explained in Sect. 13.2] for scenarios 



and in which the BBN constraints are respected, the present 
velocity of gravitinos emitted in ri decays is between 0.004 
and 0.01 km/s which is comparable to that of a thermal relic 
warm dark matter species with a mass between 1 and 5 keV; 
cf. Fig. 14 of Ref. [2^ . 

■^^ For the case of broken R-parity, see, e.g., [38ll39ll40] . 



with late-decaying G's. Indeed, these BBN constraints dis- 
favor the x\ NLSP for mg > 100 MeV [^^51[5TW7j . For 
the charged slepton NLSP case, the BBN constraints as- 
sociated with electromagnetic/hadronic energy injection 
have also been considered and found to be much weaker 
but still significant in much of the parameter space |223i 
[3T41129811671I2 07]. This can be seen for a [1 ~ [r NLSP 
with (|38p in Fig. [TOl where the constraints from electro- 
magnetic and hadronic energy release are shown respec- 
tively by the short-dashed (blue, labeled as D^J^''^""*') and 
long-dashed (green, labeled as •^He/D) lines and by the 
short-dash-dotted (blue, labeled as D^'^^''^™'') fines. For 
the [i NLSP within the CMSSM with Y^^""" calculated by 



micrOMEGAs 1.3.7 [209l[2T0]. the electromagnetic ^He/D 
and the hadronic 0^1^''^°"'^ constraints are also shown re- 
spectively by the thin gray (pink) and the short-dash- 
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TTijj = mo, tan j3 = 10, Aq = 0, /j. > 




500 1000 2000 
Tni/2 [GeV] 



5000 



Fig. 11. CMSSM regions with n^h^ G ^2^^, for Tr = 10^ 10^ 
and 10^ GeV indicated respectively by the light, medium, and 
dark shaded (green) bands in the (mi/2, w2o) planes for tan f3 = 
10, Aq = 0, n > 0, and — mo. The Tr values are associated 
with the Tr definition given in (|30|l . The regions excluded by 
the chargino and Higgs mass bounds and the line indicating 
m-o — rriri are identical to the ones shown in Fig. |3] In the 
dark gray region, the gravitino is not the LSP. The dotted lines 
show contours of the NLSP lifetime. The region to the left 
of the long-dash-dotted (red) line and to the left of the thin 
gray (pink) line is disfavored by the observationally inferred 
abundances of primordial ^Li and ^He/D [61]. The effect 
of hadronic energy injection on primordial D [298j disfavors the 
Ti NLSP region above the short-dash-dotted (blue) lines. The 
Xi NLSP region is disfavored by BBN constraints from energy 
injection |297 ,223,314,,67). On the solid vertical line (violet) 
rUg = 2.5 TeV. From [317] . 



dotted (blue) lines in Fig. \Up\ In the vi NLSP case, 
there are basically no electromagnetic constraints while 
the hadronic ones are similar to those in the charged slep- 
ton NLSP case; see Refs. [223l[224l[63l[2T9] for details. In 
the ti NLSP case, there can be electromagnetic and ha- 
dronic energy injection. Indeed, the (average) hadronic 
energy emitted in a single ti decay ehad can be relatively 
large. However, since BBN constraints apply often to com- 
binations such as eha.dYi^'^^, this can be compensated to 
some extend by a relatively small Y~°'^ , which results from 
efficient annihilation due to the strong coupling of the col- 
ored h [222ll212il308j . In fact, a ti NLSP can experience 
color confinement due to its color charge and thus allows 
for intriguing non-trivial scenarios j307ii2T^213j . 



The electromagnetic constraints shown in Figs. \W\ and [TT] 
apply to the ti NLSP and are obtained from Fig. 42 of Ref. [61] 
for a 'visible' electromagnetic energy of iJvis = eom ~ 0.3 Er_of 
the T energy Er — {rrif^ — m~-|-m^)/2m?i released in ri — » Gr. 



An additional constraint on electromagnetic energy in- 
jection can be inferred from the observed Planck spectrum 
of the CMB 321,,322J. This CMB constraint is not shown. 
According to the results of |322] for the case of a tr NLSP, 
this limit is everywhere less severe than the severe electro- 
magnetic limit Dl'^ given by the short-dashed (blue) line 
in Fig. M 

As already mentioned in the Introduction, also the 
mere presence of a long-lived negatively charged particle 
X~ — such as a long-lived /j~ — can lead to bound states 
that catalyze BBN reactions and can thereby be associ- 
ated with BBN constraints. In fact, the possibility that 
BBN can be affected by bound-state formation of an X~ 
with primordial nuclei had already been realized almost 
twenty years ago [323.324. 325] . It was however only less 
than three years ago [64 when it was realized that bound- 
state formation of with ^He can lead to a substantial 
production of primordial ^Li via the CBBN reaction 



{^IieX-) + D ^ ^U + X- 



(40) 



Since then, there has been a considerable effort to refine 
various aspects of CBBN and to understand its implica- 
tions in the framework of specific models [551,'66U671 I3091 
[207l[68l[6m [70l[7T][72l[73l[^ ,63,31 
I327j . In particular, it has been pointed out in 78J that 
there is also the possibility of efficient ^Be production via 
a radiative fusion of ^He and (*HeX~) leading to {^BeX~), 
which can capture a neutron resonantlvF^ 



^He + {'^ReX' 
CBeX-) + 



(^BeX-) +7 
^Bc + X- . 



(41) 
(42) 



The efficiency of CBBN of both ^Li and ^Be depends 
strongly on the abundance of X~ at the relevant times, 
which is given by Y-'^'^^ and rr for X^ = , i.e., in the slep- 

ton NLSP case. Observationally inferred limits on the pri- 
mordial abundances of both ^Li and ^Be can thus be used 
to extract t^^ -dependent upper limits on IC^ec^ Indeed, 

from the limits ^Li/Hlobs < 10-"-10-i° (cf. [59113281 
[73) and 9Be/H|obs < 2.1 x 10"" (79], the T,^-dependent 

bounds on Y-'^^^ shown in Fig. 5 of Ref. [79| have been 
'i 

obtained. Confronting ([38| with those bounds, one finds 
the ^Be and ^Li constraints shown respectively by the 
solid and by the long-dash-dotted (red) lines in Fig. [TU] 
The long-dash-dotted (red) line in Fig. [Tl] shows the con- 
straint associated with ^Li/Hjobs ^ 2 x 10~^^ [IS as ob- 
tained from the CBBN treatment of |74] with Y- from 
micrOMEGAs. 

For a typical yield ([38]) , the CBBN constraints associ- 
ated with ^Li and ^Be imply t64u68] l69 [ l7 H I74 p 79j 



rr < 5 X 10^ s. 



(43) 



The large ^Be-production cross section reported and used 
in Refs. [781179] has very recently been questioned by Ref. [327], 
in which a study based on a four-body model is announced as 
work in progress to clarify the efficiency of ®Be production. 
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While numerous other CBBN reactions can also affect the 
abundances of ^Li, ®Be, and other primordial elements [67t 
[69.. 73, 75, 78, 79, 327], the approximate tj^ bound is rela- 
tively robust. In particular, the possibility of allowed is- 
lands in the parameter region with large Y^^^"^ /large Ti_^ — 

which was advocated to remain viable in Ref. |73j — does 
not exist 79]. The finding of relaxed limits at long 

lifetimes tj^ in Ref. |73| did rely on the presence of [pl^)- 
bound states and on a claimed significant reprocessing of 
^Li by nuclear reactions such as (p ? f ) -I- ^Li ~* + ''He + 
^He. However, it is clarified in Ref. [79] that the pres- 
ence of (p Zj~)-bound states cannot relax the Y^'^ hmits at 
long lifetimes tj^ in any substantial way because of charge 
exchange reactions such as {pli) + *Hc {'^Hel^) -f p 
and ipli ) + ^Li (^Li/j~) -I- p that can be very effi- 
cient |79j . This finding has recently been confirmed in 
quantum three-body calculations [32 7j . 

As can be seen in Fig. [TOl the cosmological constraints 
provide an upper bound on to^ once mf^ is measured. 
This bound implies upper bounds on the SUSY breaking 
scale, on fl^^ , and — as can be seen in Fig. [9] — on Tr. 
Figure [11] shows that the cosmological constraints imply 
not only an upper limit on [207] but also a lower limit 
on mi/2 [67ll207| . Indeed, mg-dependent limits on the re- 
heating temperature, 
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Fig. 12. Upper limits on the mass ratio c = mg/mj-^ imposed 

by ^Y^^ - ■^dmfe^ < 0.126 and r,-^ < 5 x 10^ s (lO'^s) are 
shown as a function of m^-^ by the solid (dashed) lines for values 
of T]^^'' ranging from 10** GeV up to 3 x 10^ GeV. These values 
refer to the Tr definition given in (|30p . From [316] . 



Tr < 4.9 x 10 



^ GeV 



.10 GeV/ 
and on the gaugino mass parameter. 



/ \2/5 
"^V^ ^0-9 TeV (^) . 



(44) 



(45) 



have been derived within the CMSSM [Hi from the limit 
([i5|) . i.e., under the assumption that Yj is described 



approximately by ([55]) While the Tr bound can be re- 
strictive for models of inflation and of baryogenesis, the 
mi/2 bound can have implications for SUSY searches at 
the LHC. Depending on mg, ([^5]) implies sparticle masses 
which can be associated with a mass range that will be 
difficult to probe at the LHC. This is illustrated by the 
vertical (violet) line in Fig^ 11] which indicates the gluino 
mass rUg = 2.5 TeV [3T7] 



Here one should emphasize that the Tr limit re- 
lies on the CMSSM-specific minimal splitting [74] < 
0.21m^y2- Iri a less restrictive model, higher Tr values can 
be viable for a smaller splitting between the masses of the 
gluino and the h NLSP 316 . This is shown in Fig. [H] 
(from [316j ). where the maximum reheating temperature 
< f2ri„, and rr < 5 x 10^ GeV 



^^max 



imposed by 



(10** GeV) is shown by the solid (dashed) lines in the plane 



spanned by m^^ 



and the ratio c — 



at the weak 
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Similar limits have also been discussed in models where the 



ratio is bounded from below 

Note that the mass of the lighter stop is mj;_^ 
the considered ri NLSP region with mh > 114.4 GeV 



0.7mg in 



scale. Note that the shown T^'^^ values are given by [316] 
2.37 X 10^ GeV 



Tr < 



0.1 



(t^)'( 



mr 



. 100 GeV 



T^max 

R, 



(46) 



which has been derived in a very conservative way taking 
into account only the SUSY QCD contribution to 

The BBN constraints shown in Figs. [TO] [lU and [12] 
and dU), (dU), and (gS]) rely on typical Yj^"" values such 

as the ones given by ([SH]). However, F~^'='= < 4 x 10" can 
be found even within the CMSSM as is shown in Fig. [5] 
Indeed, for a ri NLSP with a sizable left-right mixing, 
exceptional values oiY^f:'^ < 10~^^ can occur |310p311j for 

which even the restrictive CBBN bounds associated with 
6Li/H|obs < 10-1° and ^Be/Hlobs < 2.1 x lO'^^ (cf. Fig. 5 
of [79 ) can be evaded. In these exceptional cases, the tj^ 
limit (|43)) does not exist and the Tr limit is governed by 
dm only (cf. Fig. [9]) since i?-"^^ is negligible. 



< n. 



4.3 Experimental Prospects 

Because of its extremely weak couplings, gravitino dark 
matter is inaccessible to direct and indirect searches if R- 
parity is conserved(3 Also the direct production of grav- 



^ For broken R-parity, the G LSP is unstable and its decay 
products can lead to signals in indirect searches |39ll329ir330l 
I33ni332ll333ll334l. 
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itinos with > 0.1 keV at colliders is extremely sup- 
pressed. Instead, one expects a large sample of (quasi-) 
stable NLSPs if the NLSP belongs to the MSSM spec- 
trum and if its mass is within the kinematical reach. 

The Xi E^rid the t'l NLSP cases will be associated with 
an excess of missing transverse energy and it might be- 
come a major challenge to distinguish such cases from the 
Xi LSP case; cf. |218[|219| . Indeed, the phenomenological 
prospects will be much more promising in scenarios with 
the /i NLSP or the ti NLSP. These particles can appear 
respectively as (quasi-) stable electrically charged elemen- 
tary particles or as stop hadrons in the collider detectors. 
Because of possibly non-trivial hadronization properties 
of the ti NLSP, we focus on the simpler ti NLSP case in 
the remainder of this section. 

In the Ti NLSP case, each heavier superpartner pro- 
duced will cascade down to the ti which will appear as 
a (quasi-) stable particle in the detector. Such a heavy 
charged particle would penetrate the collider detector in 
a way similar to muons [335.,33 6,„337^. If the produced 
staus are slow, the associated highly ionizing tracks and 
time-of-flight measurements will allow one to distinguish 
the Ti from a muon 335,336,337,338,. With measure- 
ments of the Ti velocity /Jy^ = w^^/c and the slepton 
momentum p^-^ = \pr-^ \ , ^ can be determined: mf ^ = 
p?i(l - /?- )^/V/3?i ISSH]- For the upcoming LHC exper- 
iments, studies of hypothetical scenarios with long-lived 
charged particles are actively pursued; see, e.g., 339,340, 
134111342] . For example, it has been found that one should 
be able to measure the mass m^^ of a (quasi-) stable ti 
quite accurately 339,340 13 

If some of the staus decay already in the collider de- 
tectors, the statistical method proposed in 338] could al- 
low one to measure the ti lifetime. With ([55]) and the 
measured value of m?j, one will then be able to deter- 
mine also the gravitino mass tjIq and thereby the scale 
of SUSY breaking. As a test of our understanding of the 
early Universe, it will also be interesting to confront the 
experimentally determined (m^,m^-^) point with the cos- 
mological constraints in Fig. [TOl 

Ways to stop and collect charged long-lived particles 
for an analysis of their decays have also been proposed [3441 
[3151:346,347,348^349 . It was found that up to ©(lO^-lO'') 
and O(IO^-IO^) ti's could be trapped per year at the LHC 
and the ILC, respectively, by placing 1-10 kt of massive 
additional material around existing or planned collider de- 
tectors |3451l346j . A measurement of r^^ can then be used 
to determine as already described above. If can 
be determined independently from the kinematics of the 
2-body decay ri Gt, 



= J + ml — 2m , (47) 



^® (Quasi-) stable n's could also be pair-produced in interac- 
tions of cosmic neutrinos in the Earth matter and be detected 
in a neutrino telescope such as IceCube [343j . 



the lifetime Tf ^ can allow for a measurement of the Planck 
scale [550113511134811352) 

An agreement with (|26p . which is inferred from Newton's 
constant [2] Gn = 6.709 x 10~3^GeV"^, would then pro- 
vide a strong experimental hint for the existence of su- 
pergravity in nature |350j . In fact, this agreement would 
be a striking signature of the gravitino LSP. Unfortu- 
nately, the required kinematical determination of ap- 
pears to be feasible only for |3511l3481l352j TOg/mfj ^0.1 
which seems to be disfavored in most of the SUSY pa- 
rameter space according to our present understanding of 
the cosmological constraints (see Fig. [TU| F1 Accordingly, 
alternative methods such as the ones proposed in jl90[ 
1353) could become essential to identify the gravitino as 
the LSP. In the special regions with exceptionally small 
yjiec < 20~i5 j3T0 ,311J, however, the cosmological con- 

straints could be evaded even within a standard cosmo- 
logical history and mg can still be sufficiently close to 
rn^i so that the kinematical mg determination can still 
be viable. Here, also the spin-3/2 character of the grav- 
itino can become relevant so that it could be probed in 
principle by analyzing the decays ti Gtj |350| . 

With additional experimental insights into the masses 
of the gluino and of the neutralinos (and thereby into 
-^1.2,3 and their possibly universal value toi/2 at Mqut), 
the determination of (or t^:^) can allow one to probe 
also the upper limit on the reheating temperature T^^^ at 
colliders 293i r2Mll231ll222l[M8l[23l [7Tl [3T51[3T6] : cf. Figs.H 
[TTl and [T21 This possibility results from the extremely 
weak gravitino couplings, the associated Tr dependence 
of and the limit < i7dm- Any T^"^'' value in- 

ferred from collider experiments will however depend cru- 
cially on assumptions on the cosmological history and on 
the evolution of physical parameters as discussed already 
for bounds in Xi LSP scenarios in Sect. 13.21 It would 
still be very insightful to compare, e.g., the minimmn tem- 
perature required by thermal leptogenesis in its simplest 
setting, T > 10^ GeV, with the T^'^'' value obtained under 
the assumptions of a standard cosmological history and of 
couplings that behave at high temperatures as described 
by the rcnormalization group equations in the MSSM. 

5 Axino Dark Matter 

The axino a j354p355)l356p357) appears as the spin- 1/2 su- 
perpartner of the axion once the MSSM is extended with 
the PQ mechanism [2T|[22] in order to solve the strong CP 
problem. Depending on the model and on the SUSY break- 
ing scheme, the axino mass mj; can range between the eV 

Note that the cosmological constraints described in 
Sect. IT2I assume a standard thermal history. In fact, entropy 
production after NLSP decoupling and before BBN can weaken 
the BBN constraints significantly [31311207] . 
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and the GeV scale j356 p 358 p 359 p 360 p 36i p 362j . The axino 
is a singlet with respect to the gauge groups of the Stan- 
dard Model. It interacts extremely weakly since its cou- 
plings are suppressed by the PQ scale /a ^ 6 x 10^ GeV [21 
[TT?[|^I^ . and thus can be classified as an EWIP. The 
detailed form of the axino interactions depends on the ax- 
ion model under consideration; cf. Sect. [2l We focus on 
hadronic (or KSVZ) axion models |,102..103J in a SUSY 
setting, in which the axino couples to the MSSM particles 
only indirectly through loops of additional heavy KSVZ 
(s)quarks. Considering a LSP scenarios in which the LOSP 
is the NLSP, m^j^ < m^o is again viable as also the alter- 
native Xij i^i: and ti NLSP cases. 

Before proceeding, it should be stressed that the bo- 
sonic partners of the axino, the axion and the saxion, can 
have important implications for cosmology: (i) Those as- 
sociated with the axion have already been discussed in 
Sect. [2] and most important for this section is the possibly 
significant axion contribution to ^7dm which can tighten 
the constraints from J7j^ < i7dm discussed below, (ii) Late 
decays of the saxion can lead to significant entropy pro- 
duction 364, 365, 366,367^ and can thereby affect the cos- 
mological constraints |368j . In this review, however, a stan- 
dard thermal history is assumed which implies that saxion 
effects are negligible. 
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Fig. 13. Upper limits on the reheating temperature Tr in 
the a LSP case for fa/N = lO" GeV. On (above) the gray 
band, n'^^h^ € 0.113l{j;gl| (J^J^/i^ > 0.129). Thermally pro- 
duced axinos can be classified as hot, warm, and cold dark 
matter [371] as indicated. From |369| . 



5.1 Primordial Origin 

Because of their extremely weak interactions, the tem- 
perature Tf at which axinos decouple from the thermal 
plasma in the early Universe is very high. For example, 
an axino decoupling temperature of Tf « 10^ GcV is ob- 
tained for fa = 10" GeV [355,369J. For Tr > Tf, axinos 
were in thermal equilibrium before decoupling as a rela- 
tivistic species so that |359 ll370l[37T1 [369 ] 

where g*s{Ti) — 228.75 is assumed. 

For Tr < Tf, axinos have never been in thermal equi- 
librium with the primordial plasma but can be generated 
efficiently in scattering processes of particles that are in 
thermal equilibrium with in the hot MSSM plasma [370[ 
[37111369113721 . Within SUSY QCD, the associated ther- 
mally produced (TP) axino density can be calculated in 
a consistent gauge-invariant treatment that requires weak 
couplings (5s < 1) |369j : 

r.TP,2 . ^ 6/^ M f 1-108 \ /10"GeV\' 

><fe)(io^) 

with the axion-model-dependent color anomaly N of the 
PQ symmetry. The thermally produced axinos do not af- 
fect the thermal evolution of the LOSP (or NLSP) which 
decays after its decoupling into the a LSP. Taking into 



account the non-thermally produced (NTP) density from 
NLSP decays [5751I57T] 

nr^'h' = m~a Y^ll^ s(To)ftVPc , (51) 
the guaranteed axino density 

= ^„thco„/TP ^ ^NTP (52) 

In Fig.dKfrom [3MI) the (m^, Tr) region with 0.097 < 
QJP < 0.129 for fa/N'^ 10" GeV is shown by the gray 
band. Note that ([50]) shows a different dependence on the 
LSP mass than the corresponding expression in the G LSP 
case (p4l) . Accordingly, one finds the different ttilsp depen- 
dence of the Tr limits inferred from illY^ < i^dm- Since 

a/G 

thermally produced axinos are generated in kinetic equi- 
librium with the primordial plasma, they have a thermal 
spectrum which allows for the m^-dependent classifica- 
tion into cold, warm, and hot dark matter |371j shown in 
Fig. [131 As can be seen, the Tr limit does not exist for 
rria ^0.2 keV because of the equality of a production and 
a disappearance rates for T > T{ ^ 10^ GeV. With a ther- 
mal relic density (|49p in this regime, there will be a limit 
on rria depending on the constraints inferred from studies 
of warm/hot dark matter [7]. 

The non-thermally produced axino density S^-"^^ dif- 
fers from the corresponding expression in the G LSP case 
P6p only by the obvious difference in m-^^. In particu- 
lar, for given fi~^ , the bound f^~^^ < f^dm — ^u'^^"^^ as 
obtained with (|38|) implies limits on rria and m^^ which 

Axino production in inflaton decays is not considered. 



Frank Daniel Steffen: Dark Matter Candidates 



23 



can be read off directly from Fig. [10] after the replace- 
ment itIq — > rria- Note, however, that the r^^j^ contours 
and the cosmological constraints are different in the ax- 
ino LSP case. For the ri NLSP, the following lifetime was 
estimated [190] 



where = 2E^/mf-^ is the scaled photon energy, 9 the 
opening angle between the directions of 7 and r, 



l-Ai 



25 sr 



dx-yd cos 9 



dcost 



(55) 



^ 1^ 100 GeV Y /JCaYY '^f 100 GeV Y ^^3^ 



TTif, A 10" GeV y 



B 



with the KSVZ-model dependence expressed by CaVY — 
0(1) and the uncertainty of the estimate absorbed into 
^ ~ 0{1). One thus finds a ti lifetime in the a LSP case 
that cannot be as large as the one in the G LSP case (p9|) . 
Accordingly, the BBN constraints are much weaker for 
the a LSP. For discussions of a LSP scenarios with the 
Xi NLSP, see [87811371 Il374j . For both the n NLSP and 
the Xi NLSP, it has been shown that non-thermally pro- 
duced axinos with < 10 GeV would be warm/hot dark 
matter ^320) . 



5.2 Experimental Prospects 

Being an EWIP, the axino LSP is inaccessible in direct and 
indirect dark matter searches if R-parity is conserved. Also 
direct a production at colliders is strongly suppressed. 
Nevertheless, (quasi-) stable ri's could appear in collider 
detectors (and neutrino telescopes [343j ) as a possible sig- 
nature of the a LSP. However, since the Mp measurement 
at colliders [350j , which would have been a decisive test of 
the G LSP, seems cosmologically disfavored in most of the 
parameter space of R-parity conserving SUSY models, it 
may very well be a challenge to distinguish between the a 
LSP and the G LSP. 



the respective integrated 3-body decay rate with the cuts 
x^ > x^y"' and cos 9 < 1 — Xg"' , and Ai = m,f / m?^ . Note 
that ([54]) is independent of the 2-body decay, of the total 
NLSP decay rate, and of the PQ/Planck scale. The figure 
shows (|5^ for the axino LSP (i — a) with m~/m~^ <C 1 

(upper panel) and the gravitino LSP {i = G) with = 
10 MeV (lower panel), where nif-^ ~ 100 GeV, nij^ — 

110 GeV, and x™* = 4"' = 0.1. In the G LSP case, the 
events are peaked only in the region where photons are soft 
and emitted with a small opening angle with respect to 
the tau {9 ~ 0). In contrast, in the a LSP case, the events 
are also peaked in the region where the photon energy is 
large and the photon and the tau are emitted back-to-back 
(0 ~ tt). Thus, if the observed number of events peaks in 
both regions, this can be evidence against the gravitino 
LSP and a hint towards the axino LSP [l90]FI 

To be specific, with 10'* analyzed stau NLSP decays, 
one expects about 165±13 (stat.) events for the a LSP and 
about lOOilO (stat.) events for the G LSP 190J, which 
will be distributed over the respective (x-y, cos^)-planes 
shown in Fig. 1141 In particular, in the region of Xj > 0.8 
and cos 9 < —0.3, we expect about 28% of the 165±13 
(stat.) events in the a LSP case and about 1% of the 
lOOiblO (stat.) events in the G LSP case. These numbers 
illustrate that 0(10^) of analyzed stau NLSP decays could 
be sufficient for the distinction based on the differential 
distributions. To establish the feasibility of this distinc- 
tion, dedicated studies including details of the detectors 
and the additional massive material will be crucial |348j . 



For 



100 GeV and = 110 GeV, for exam- 



5.3 Probing the Peccei-Quinn Scale fa and 

If a is the LSP and ti the NLSP, the analysis of the 2- 
body decay ti to, will allow us to probe the PQ scale 
fa and the axino mass m^. In fact, the measurement of 
T^:^ (l53l) with methods described in Sect. 14.31 leads to the 



pie, the Ti lifetime in the a LSP scenario (|53[) can range 
from 0(10-* s) for fa ^ 5 x 10^ GeV to 0(10 h) for 
fa = 5^ 10*'^ GeV. In the G LSP case, the corresponding 
lifetime (j39|) can vary over an even wider range, e.g., from 
6 X 10^^ s for TOg — 1 keV to 15 years for mg = 50 GeV. 
Thus, both a very short lifetime, Tr^ < ms, and a very long 

one, Tr, > days, wiU point to the G LSP. On the other foUo'whig estimate of the Peccei-Quinn scdT/a [l90]: 
hand, if the LSP mass cannot be measured kinematically 
and if t^-^ = 0(0.01 s)-O(10 h), the stau lifetime alone 
will not allow us to distinguish between the a LSP and 
the G LSP. 

The situation is considerably improved when one con- 
siders the 3-body decays ri r^a/G. From the corre- 
sponding differential rates [190] . one obtains the differen- 
tial distributions of the visible decay products. These are 
illustrated in Fig. [14] (from [190j ) in terms of 



fl^ecl^Y{l^''GeY)\l~ 

X ( "^^1 ^ ( "^B \ 

V 100 GeV/ V 100 GeV/ 



25 s 



(56) 



which can be confronted with limits from axion stud- 
ies; cf. Sect. [2] Indeed, we expect that m^^ and will 



1 



dx~fd COS 9 



(54) 



^® There is a caveat: If mg < ma < m^^ and -r(ri — > aX) 3> 
_r(ri — + GX), one would still find the distribution shown in 
the upper panel of Fig. 1141 The axino would then eventually 
decay into the gravitino LSP and the axion. 
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Fig. 14. The normalized diflterential distributions (|54|l of the 
visible decay products in the decays ri t + 'y + a/G for the 
cases of the a LSP (upper panel) and the G LSP (lower panel) 



for m? 



100 GeV, 
= lOMeV 



X? 



m; 



= 110 GeV, ml/mi < 1 



The cut parameters are set to x'. 



cut 

7 



and mg 

ig"' — 0.1. The contour lines represent the values 0.2, 0.4, 0.6, 
0.8, and 1.0, where the darker shading implies a higher number 
of events. From 11901. 



already be known from other processes when the ti NLSP 
decays are analyzed; cf. Sect. 14.31 The dependence on rria 
is negligible for rria/mf-^ < 0.1. For larger values of m^, 
the Ti NLSP decays can be used to determine toj from the 
kinematics of the 2-body decay, i.e., from a measurement 
of the energy of the emitted tau Er , 



(57) 



with an error governed by the experimental uncertainties 
on and Er- As is evident from (|50p and ((5T|) . the 



determination of both the PQ scale fa and the axino mass 
nia is crucial for insights into the cosmological relevance 
of the axino LSP. 

In principle, the determination of both fa and tos at 
colliders would allow one to probe an upper limit on the 
reheating temperature T^'^^ at colliders for ma > 0.2 keV 
as can be seen in Fig. [131 see also j315j . Indeed, Ref. |315j 
has outlined theoretical ways to probe Tr values as low as 
100 GeV with fi-"^^ taken into account but without ad- 
dressing the difficulties to determine an axino mass < 
10 GeV experimentally. In analogy to the G LSP case, 
the theoretical possibility results from the extremely weak 
axino couplings, the associated Tr dependence of f^~^ , 
and the limit f]~^ < i?dm- However, the high Tr values 
at which the gauge-invariant result for f2~^ is reliable, 
Tr > 10^ GeV, are associated with rria <C 1 GeV while the 
kinematical determination ()57p appears to be feasible only 
for rria > O.lrrif, > 10 GeV, as in the G LSP case [35T| 
I348[l352j . Since also t^-^ is practically TOj-independent for 
'^o ^ 10 GeV < O.lm^j, it seems unfortunately to be im- 
possible to probe those Tr values in the axino LSP case 
at colliders. 



6 Conclusion 

The existence of dark matter provides strong evidence for 
physics beyond the Standard Model. Extending the Stan- 
dard Model with PQ symmetry and/or SUSY, the axion 
and/or an electrically neutral and color neutral LSP ap- 
pear as promising dark matter candidates. The axion is 
well motivated by the PQ solution to the strong CP prob- 
lem. With and without SUSY being realized in nature, the 
axion can exist and can contribute significantly to fidm- 
In SUSY extensions of the Standard Model, the lightest 
neutralino Xij the gravitino G, or the axino a can be the 
LSP and as such explain the non-baryonic dark matter 
in our Universe. The neutralino Xi is already part of the 
MSSM which provides a solution of the hierarchy problem 
and allows for gauge coupling unification. Being the su- 
perpartner of the graviton and the gauge field associated 
with supergravity, the gravitino G is equally well moti- 
vated with a mass mg that reflects the SUSY breaking 
scale. As the superpartner of the axion, also the axino a 
appears naturally once the strong CP problem is solved 
with the PQ mechanism in a SUSY setting. 

While mass values and interactions can be very differ- 
ent for the a, the Xi LSP, the G LSP, and the a LSP, I 
have illustrated for each of these dark matter candidates 
that there are natural regions in the associated parame- 
ter space in which fia/LSP — ^dm- In the axion case, the 
f2a = f^dm region can be subject to constraints from limits 
on axionic isocurvature and non-Gaussian perturbations 
inferred from the CMB anisotropics. In the SUSY case, re- 
gions with /2lsp = ■!^dm are limited most importantly by 
bounds from electroweak precision observables, B-physics 
observables, Higgs and sparticle searches at LEP, and by 
BBN constraints. The constraints from fidm and BBN also 
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Table 2. Dark matter candidates, their identity, and key properties. With the listed production mechanisms, Qx = Odm 
is possible for each candidate X. The respective production mechanism leads typically to a cold, warm, or hot dark matter 
component as indicated. Quantities marked with '(?)' seem to be unaccessible in large parts of the parameter space in light of 
the current understanding of experimental feasibility and/or of cosmological constraints within a standard thermal history. 



candidate 


identity 


mass 


interactions 


production 


constraints 


experiments 


a 


axion 


< 0.01 eV 


iP/fa)" 


misalign, mech. 


^ cold 


direct searches with 




(spin 0) 




extremely weak 






microwave cavities 




N.-Goldst. boson 




/a > 6 X 10^ GcV 






^ ^ ^a, /a, 5a77 




PQ symm. break. 








CMB 






lightest neutralino 


C(iOO OeV ) 


g. g'. Vi 


therm, relic 


< — cold 


indirect searches 




(spin 1/2) 




weak 


G decay 


^ warm/hot 


direct searches 




mixture of 




A/^ ^ 100 GcV 






collider searches 




B, W, Hi H'^ 








BBN 


^ n^x? ' X? coupl. 


G LSP 


gravitino 


eV-TeV 


(p/Mp)" 


therm, prod. 


<— cold 


Ti prod, at colliders 




(spin 3/2) 




extremely weak 


NLSP decay 


^ warm 


+ Ti collection 




superpartner 




Mp = 2.4 X 10^^ GcV 






+ ri decay analysis 




of the graviton 








BBN 


^ mg, Mp (?), Tr 


a LSP 


axino 


eV-GeV 




therm, prod. 


^ cold/ warm 


Ti prod, at colliders 




(spin 1/2) 




extremely weak 


NLSP decay 


^ warm/hot 


+ ri collection 




superpartner 




/a > 6 X 10^ GeV 






+ ri decay analysis 




of the axion 








BBN 


^ ma (?), fa, Tr (?) 



imply restrictive upper limits on the reheating tempera- 
ture after inflation Tr which can be relevant for models of 
inflation and of baryogenesis. 

Most promising are the experimental prospects in the 
case of the x? LSP. Being a WIMP, the x? LSP should 
be accessible in direct and indirect dark matter searches. 
Indeed, first hints might have already been found in the 
EGRET data 259 \. With ongoing indirect searches, the 
increasing sensitivity of direct searches, and the advent of 
the LHC at which xi dark matter could be produced, we 
will be able to test whether these hints are indeed the first 
evidence for the existence of SUSY dark matter. While an 
excess in missing transverse energy could provide a first 
hint for SUSY at the LHC already within the next three 
years, the identification of the Xi being the LSP will re- 
quire the reconstruction of the SUSY model realized in 
nature. If superparticles are within the kinematical reach, 
precision studies at the ILC will be crucial for this en- 
deavor. 

Also an axion discovery is conceivable to occur in the 
near future. Indeed, a current upgrade of the microwave 
cavity experiment AD MX aims at a significant experimen- 
tal exploration of a "natural" part of the Qa — ^dm re- 
gion within the upcoming years. Complementary to that, 
axion searches with helioscopes such as CAST and those 
that are performed purely in the laboratory will extend 
their search ranges. However, a signal in those searches 
would support the existence of axions as hot/ warm dark 
matter which can only provide a minor fraction of fidm- 
This would leave room for a significant f2hsp contribution 

to J7dm- 



In G/a LSP scenarios with conserved R-parity, no dark 
matter signal should appear in direct or indirect searches. 
However, since an electrically charged lightest Standard 
Model superpartner such as the ti is viable in G/a LSP 
scenarios, (quasi-) stable ri's might occur as muon-like 
particles instead of an excess in missing transverse energy. 
Indeed, an excess of (quasi-) stable ri's could appear as 
an alternative first hint for SUSY at the LHC in the next 
three years. Because of the severe limits on the abundance 
of stable charged particles [2j , one would then expect that 
the Ti is the NLSP that decays eventually into the G/a 
LSP or that R-parity is broken. A distinction between 
those scenarios will require the analysis of ti decays. For 
this challenge, the ILC with its tunable beam energy seems 
crucial [M57346 SSTl lMSllMlISp . 

Table [2] presents an overview of the dark matter can- 
didates discussed in this review. The axion and each LSP 
candidate — the lightest neutralino Xi, the gravitino G, or 
the axino a — could provide fidm and could be produced 
and identified experimentally in the near future. 
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